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We calculate finite-size effects of the Gaussian model in a L x L"'"^ box geom- 
etry with free boundary conditions in one direction and periodic boundary 
conditions in d — 1 directions for 2 < d < 4. We also consider film geometry 
{L — i> oo). Finite-size scaling is found to be valid for d < 3 and li > 3 but log- 
arithmic deviations from finite-size scaling are found for the free energy and 
energy density at the Gaussian upper borderline dimension d* = 3. The log- 
arithms are related to the vanishing critical exponent 1 — a — u = (d — 3)/2of 
the Gaussian surface energy density. The latter has a cusp-like singularity in 
d > 3 dimensions. We show that these properties are the origin of nonscaling 
finite-size effects in the mean spherical model with free boundary conditions 
in d > 3 dimensions. At bulk Tc in d = 3 dimensions we find an unexpected 
non-logarithmic violation of finite-size scaling for the susceptibility x ~ -^^ 
of the mean spherical model in film geometry whereas only a logarithmic 
deviation x ~ L^lnL exists for box geometry. The result for film geometry 
is explained by the existence of the lower borderline dimension di = 3, as 
implied by the Mermin- Wagner theorem, that coincides with the Gaussian 
upper borderline dimension d* = 3. For 3 < < 4 we find a power-law 
violation of scaling x ~ L'^^^ at bulk for box geometry and a nonscal- 
ing temperature dependence Xsurface ~ C,'^ of the surface susceptibility above 
Tc. For 2 < d < 3 dimensions we show the validity of universal finite-size 
scaling for the susceptibility of the mean spherical model with free boundary 
conditions for both box and film geometry and calculate the corresponding 
universal scaling functions for T > Tc. 
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I. Introduction and summary 



Finite-size effects near phase transitions and the concept of finite-size scahng 
near critical points have been the subject of many studies over the past 
decades |I|, ^ ^, ^. Consider, for example, the susceptibility xi^jL) of a 
ferromagnetic system at the reduced temperature t = (T — Tc)/Tc > near 
the bulk critical temperature in a cubic geometry with a linear size L 
below the upper critical dimension = 4. The property of finite-size scaling 
means that, for sufficiently large L and small t, x has the asymptotic form 

x{t,L) = x(t,oo)/^(L/0 (1) 

where x(^) oo) = AyJ:~'^ is the bulk susceptibilty and ^ = ^o^"'^ is the bulk 
correlation length. An appealing feature of finite-size scaling is universality 
which means that all nonuniversal parameters of the confined system can 
be absorbed entirely in the bulk amplitude and in the bulk correlation 
length ^, thus finite-size scaling functions such as /^(a^) are expected to be 
independent of nonuniversal details (such as the lattice structure, the lattice 
spacing and the magnitude of coupling constants). This implies that the 
amplitude of the small-x behavior f(x) = B^x"^'/" for T — at fixed L 
is also universal. The specific shape and the amplitude B^ of such scaling 
functions do, of course, depend on the geometry and on the kind of boundary 
conditions. A central prediction of finite-size scaling is the size dependence 
at the bulk critical temperature Tc 

X(0,L) = A^i^''^ B^L^I'^ (2) 

with the hulk critical exponent 7/z/ and the universal amplitude B^. For 
purely periodic boundary conditions and short-range interactions, universal 
finite-size scaling in the sense of Eqs. (|l]) and @) has been largely confirmed, 
except for the nonnnroersa/ exponential behavior in the region L » ^ which 
has recently been shown p, |^, |^ to depend on the lattice structure for 
lattice models and on the cutoff procedure for continuum models. 
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Of particular interest are nonperiodic boundary conditions which are relevant 
for real systems. For example, for the superfiuid transition of ^He, Dirichlet 
boundary conditions of field theories are believed to be fairly realistic p. 
For this system, however, accurate experiments have detected nonscaling 
finite-size j^, |10|, |Tl], |12[ and surface [1^ effects that dbTQ clS yet unexplainec. 



Furthermore there exist unexplained finite-size effects in the XY model with 
nonperiodic boundary conditions as detected by Monte Carlo simulations 

m 



On the theoretical side, the true conditions for the validity of universal finite- 
size scaling for systems with nonperiodic boundary conditions are not estab- 
lished. This includes the important case of free boundary conditions for 
lattice models which are believed to be asymptotically equivalent to Dirich- 
let boundary conditions of continuum models. It is known that universal 
finite-size scaling in the sense of Eq. (1) fails for the mean spherical model in 
film geometry with free boundary conditions in d = 3 and d = 4 dimensions 
0, [1^, , and similarly for the ideal Bose gas with Dirichlet boundary 



conditions for cubic and film geometry 0, |T8|, In these models the 

bulk correlation length could not be used as the only reference length and 
nonscaling finite-size effects were incorporated in nonuniversal shifts of the 



temperature variable 0, 16, f7l nsl O, EO]. Logarithmic nonscaling finite 



size effects that depend on the lattice spacing exist also in Gaussian interface 



models |21[] as well as in other models 



On the other hand, universal amplitude ratios have been found for critical 
systems contained in parallel plates with nonperiodic boundary conditions 



||22|| . Furthermore, field-theoretic renormalization-group calculations have 
apparently confirmed the validity of universal finite-size scaling within the 
(f^ field theory with Dirichlet boundary conditions: Universal finite-size am- 
plitude ratios and universal finite-size contributions to the free energy 
density and to the critical Casimir force were calculated both in the Gaus- 



sian (one-loop) approximation [U, ^ ^ as well as in two-loop order ||25| 
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Universal finite-size scaling functions have also been predicted for the specific 
heat and the superfiuid density in the presence of Dirichlet boundary condi- 
tions P, Related field-theoretic predictions have also been presented for 
surface quantities |2|. In these papers ^ |2|, ^ , how- 



ever, the method of dimensional regularization was employed which neglects 
lattice and cutoff effects. 

Recent work on finite-size effects ||^, ^, ^ ^ ^ |S3|, 0] has demon- 
strated that general renormalization-group arguments are not sufficient to 
prove the validity of universal finite-size scaling and that cutoff and lat- 
tice effects are nonnegligible for confined systems with periodic boundary 
conditions. Clearly these investigations need to be extended to the case of 
nonperiodic boundary conditions. 

The corresponding analytic calculations, at finite cutoff and at finite lattice 
spacing, become quite difficult within the (p^ theory beyond the lowest order. 
Before embarking on such an ambitious project it is of course necessary to first 
examine the lowest-order case under the simplest nontrivial conditions, i.e., 
with free (or Dirichlet) boundary conditions in only one direction. Therefore, 
as a first step, we consider the exactly solvable Gaussian model with short- 
range interaction on a simple-cubic lattice with a lattice constant a for a 
finite rectangular L x L'^~^ box geometry with free boundary conditions in 
one direction and periodic boundary conditions in d — 1 directions. Even at 
the Gaussian level, the analytic calculations at finite lattice spacing in the 
range 2 < d < 4 turn out to be nontrivial. 

For the specific heat and the susceptibility of the Gaussian model we find 
full agreement with universal finite-size scaling. With regard to the singular 
part of the free energy we find that the finite-size scaling form is indeed 
valid for d < 3 and d > 3 but logarithmic deviations from finite-size scaling 
occur at d = 3 where the critical exponent 1 — a — z/ = (c? — 3)/2 of the 
surface energy density vanishes. In order to describe the logarithmic d = 3 
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behavior it is necessary to keep the lattice spacing finite. We find that the 
same logarithmic deviations from finite-size scaling exist in the continuum 
version of the Gaussian model with Dirichlet boundary conditions provided 
that a finite cutoff is used. This implies that the method of dimensional 
regularization at infinite cutoff is not capable of correctly describing the 
d = 3 behavior of the singular part of the free energy density and of the 
energy density since it yields unphysical divergences of these quantities in 
the form of a pole term ~ (c? — 3)~^ ||35|. As discussed in Sect. III. H, the 
dimension d = 3 can be considered as an upper borderline dimension d* of 
the Gaussian model with free boundary conditions above which lattice and 
cutoff effects become nonnegligible for the surface energy density. 

For d > 3 we find that the surface energy density Usurfaceit) of the Gaussian 
model with free boundary conditions has a cusp-like singularity at bulk Tc as 
Tc is approached from above. For the lattice model at finite lattice spacing 
d the height of the cusp is 

with 



oo 



B, = ^ldy{[l + e-'y - 2 e~'y hi2y)] [e'^y hi^y)]'''] > (4) 



where loiz) is the Bessel function of order zero. The temperature dependent 
part of Usurfaceit) has a universal scaling form ~ but it vanishes at 
and is sub leading compared to the nonuniversal finite regular part, Eq. @, 
at Tc. The latter part yields a leading nonscaling contribution 2 Usurfacei^) / L 
to the total energy density. These results remain valid also for the Gaussian 
model in film geometry (L oo) with free boundary conditions. 

In a second step we analyze the exactly solvable mean spherical model with 
the same boundary conditions. Previously this model has been studied for 



film geometry at integer dimensions d = 3,4,5, ... ]T5| |T^. Here we extend 
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this analysis to continuous dimensions in the range 2 < d < A and consider 
both film and box geometry. This reveals (i = 3 as a borderline dimension 
between a universal scaling {d < 3) and a nonuniversal nonscaling [d > 3) 
regime. In this paper we calculate the nonscaling effects for 3 < < 4 as 
well as the analytic form of the universal finite-size scaling function f^{L/^), 
Eq. (|l]), of the susceptibility for 2 < c? < 3 including the amplitude B{s) of 
the scaling result, Eq. (Q) with 7/z/ = 2 , 

X(0,L) = Bis)L' , d<3 (5) 

at arbitrary shape factor s = L/L > 0. The amplitude B{s) is shown to 
diverge for (i — 3. 

The mean spherical model can be considered as a Gaussian model with a 
constraint where the constraint can be expressed in terms of the Gaussian 
energy density. Our results for the latter quantity explain the origin of loga- 
rithmic nonscaling terms in thermodynamic quantities of the mean spherical 
model ai d = 3 and of power-law violations of finite-size scaling for 3 < d < 4. 
While previous work suggested the existence of only logarithmic deviations 
from finite-size scaling in = 3 dimensions 0, |1^, |16|, |1^, |1^, ^ we 
find, quite unexpectedly, a non-logarithmic violation of the scaling predic- 
tion, Eq. (^, for the size dependence of the susceptibility at bulk Tc in c? = 3 
dimensions 

X(0,L) = Bfumd-'L^ (6) 

for film geometry whereas for box geometry, at fixed finite shape factor s = 
L/L > 0, we find the expected logarithmic deviation from scaling 

x{0,L) = B,Us)L'HL/a) . (7) 

As will be shown in detail in Sect. IV. C, the special result of Eq. (P) for film 
geometry at d = 3 is due to the simultaneous appearance of two logarithmic 
effects at d = 3 where two borderline dimensions coincide : it is a combined 
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effect of tlie logaritfimic surface term of the Gaussian model at the (upper) 
borderhne dimension d* = 3 where the exponent 1 — a — u vanishes and 
of a logarithmic finite-size term arising from the mode continuum of the 
film system just at the (lower) borderline dimension rf; = 3 at which the 
film critical temperature vanishes in accordance with the Mermin- Wagner 
theorem ||3^. Most striking is the discontinuous change of the exponent 2 
of the power law XfUm = -B(0)L^ for d < 3, Eq.(||), to 3 of the power law 
Xfiim = Bfiim d-^L^ for d = 3, Eq. (1). 

The result of Eq. (^ is not contained in the work of Barber and Fisher |1T5[| 
who calculated x for film geometry in d = 3 dimensions only for T > T(L) 
where T(L) > is some temperature that they called " quasicritical" . Our 
d = 3 result for x(t, L) covers the entire critical region T > Tc including the 
regime T > T{L). In the latter regime, the explicit form of our result is at 
variance with the simpler form of Barber and Fisher. 

For box geometry in 3 < d < 4 dimensions we find a power-law violation of 
scaling at Tc 

X(0,L) = BboAs,d)d'-''L''-' (8) 

where the amplitude Bhox{s, d) is proportional to the amplitude B^, Eq. (H), 
of the cusp of the Gaussian surface energy density. A nonscaling form is 
also found for the temperature dependence of the surface susceptibility for 
3 < d < A above Tc, 

y f = A f a^-^ f'^ ~ t-'^/id-2) (Q\ 

with ~ t^" ,v = {d — 2)^^ , whereas the scaling form, Eq. (|l]), would 
imply Xsurface ~ 0{xb ~ t^3/(d-2) ^^^q mean spherical model. Again 
the amplitude Agurface in Eq. (^) is proportional to Bd. 

For film geometry in d > 3 dimensions we find an anomalous enhancement 
of the film critical temperature T^^diL) above the bulk critical temperature 
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Tc^rf(oo). A corresponding shift was first found for c? > 4 by Barber and 



Fisher |T^. This enhancement is most naturally expressed in terms of the 
dimensionless parameter 2jPc,d{L) = 2J[kBTc^d{L)]~^ where J is the nearest- 
neighbor coupling. The result is for d > 3 

2J (oo) - /3,,, (L)] = AB, d/L - {d/Lf-^ + Oid"'^ L""/^) (10) 

with the nonuniversal amplitude Bd-, Eq. (^), and with a universal ampli- 
tude Cd > 0. Eq. (pUj) implies Tc^d{L) > Tc^rf(oo) for large L ^ d. The 
leading term ~ in Eq. ( p!OD has a nonscaling L dependence whereas the 
subleading universal term has the scaling L dependence ~ L^^'^ . 

In summary we see that both the anomalous nonscaling enhancement of 
Tc^d{L), Eq. (p!OD, and the power-law violations, Eqs. (|D and @, for c? > 3 
can be traced back to the same amplitude Bd, Eq. (^), of the nonscaling cusp 
of the Gaussian model. Thus the analysis of the Gaussian model provides 
a better understanding of the origin of the power-law nonscaling finite-size 
effects in the mean spherical model for d > 3, and, for box geometry, of the 
logarithmic deviations at the Gaussian upper borderline dimension d* = 3. 
For film geometry, however, the Gaussian logarithmic effect at d* = 3 is 
enhanced by a second logarithmic effect due to the lower borderline dimension 
di = 3 (where the film critical temperature vanishes), which then yields the 
power law Eq. (|^). 

We point out that all nonscaling effects are tied to the finite lattice constant 
a > 0, as seen explicitly in Eqs. (^ and (^ - ([T0|). We expect that similar 



effects exist in the ideal Bose gas with Dirichlet boundary conditions |]T8|, |T^, 
^ with a finite cutoff (even if a smooth cutoff is used) . These effects are not 
captured by the standard method of dimensional regularization. It remains 
to be seen whether the mechanism for nonscaling finite-size effects in the 
mean spherical model and the ideal Bose gas is an artifact restricted to these 
models or whether some of these features are of more general significance. 
This question is of particular interest below where an explanation of the 
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pronounced nonscaling finite-size effects in He remain to be a challenge for 
future research. 

In Section II we summarize the predictions implied by the finite-size scal- 
ing hypothesis. Section III contains the detailed results for the finite-size 
effects in the Gaussian lattice model with free boundary conditions and in 
the Gaussian continuum model with Dirichlet boundary conditions. In Sect. 
IV we analyze the consequences of our results for the mean spherical model 
with free boundary conditions. The derivation of our results is presented in 
several Appendices. 
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II. Finite-size scaling predictions 



In the subsequent sections we shall present exact results for the finite-size 
effects on the free energy density, energy density, specific heat and suscep- 
tibility of lattice models in a rectangular L x L*^"^ box geometry with free 
boundary conditions in the direction of size L and periodic boundary condi- 
tions in the d — 1 directions of size L. For the sake of clarity we first sum- 
marize the predictions implied by the finite-size scaling hypothesis which, 
for this geometry and these boundary conditions, have not yet been formu- 
lated explicitly in the literature. We denote the critical temperature of the 
d-dimensional bulk (L oo, L —>■ oo) system by T^^^. In the limit L oo at 
fixed L, the box becomes a film of thickness L which may have its own critical 
temperature Tc^d{L) ^ T^^^ = Tc^d(oo). In general one expects T^^d^L) < Tc^d 
but it turns out (see Sect. IV, see also Ref. |jT5|]) that for the mean spheri- 
cal model with d > 3 the film critical temperature Tc^d{L) exceeds the bulk 
critical temperature T^^d- For simplicity, in this Section, we assume a d di- 
mensional box with a finite shape factor L/L > and confine ourselves to 

First we consider the free energy density f{t,L,L) (in units of ksT) at the 
reduced temperature t = {T — Tc^d)/Tc,d > and at vanishing external field. 
It is expected that, for small t, f can be decomposed into a singular and a 



" nonsingular" part fS^, q8 



f{t,L,L) = fs{t,L,L) + fns{t,L,L) (11) 

where fns{t,L,L) has a regular t dependence. In the bulk limit the corre- 
sponding decomposition is 

m = fit, oo,oo) = fut) + m (12) 

where the regular part fo(t) = fns(t, oo, oo) can be identified unambiguously. 
For systems with short-range interactions below the upper critical dimension 
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d = 4: and for large L, L and ^ it is expected that the singular part fs{t, L, L) 
has the finite-size scaling form |l|, |3^ 

fsit.L.L) = L-^J^iL/i.L/L) (13) 

where ^(t) = ^o^"'^ is the (second- moment) correlation length of the d- 
dimensional bulk system. For a given shape factor s = L/L, the scaling 
function JF(x, s) is expected to be universal. More specifically, the singular 
and nonsingular parts of the free energy density are expected to have the 
asymptotic (small t, large L, large L) form ^, ^ 

Mt, L, L) = + 2 + <3{L/i, L/L) (14) 

and 

fns{t,L,L) = fo{t) + 2^^{t)/L (15) 

with a universal bulk amplitude and a universal surface amplitude A^urface^ 
and with a universal finite-size part Q{L/^,L/L) of the scaling function 

^(L/e, L/L) = R^iL/O" + 2 At^^f^jL/0'-' + 5(1/^, L/L) . (16) 
Eqs. (0) - (|TB|) imply that there exists the surface free energy 



L— >oo 



L 



fsurfaceit) = \im ^ f{t,L,L) - Mt) -\ = A^^. i'-" + ^l(t) (17) 



surface 



with a universal amplitude ^^^j^ce singular part. The nonsingular 

surface contribution \l/i(t) is a re^ftx/ar function of t. Nonasymptotic Wegner 
]39| corrections to scahng are neglected in Eqs. ([13|) , (p!^) and ([T6|) . The 



phenomenological finite-size scaling theory does not make specific predictions 
about the dependence on L and L of higher-order terms in Eq. (plSf). 

Eqs. ([TT|) - ([T7|) are expected to hold also for continuum models with Dirich- 
let boundary conditions in one direction, with the same universal quantities 
as for free boundary conditions of lattice systems. As noted in the Intro- 
duction, however, there exist nonuniversal exponential terms in the regime 
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L ^ ^, L ^ ^, where the lattice-dependent and cutoff-dependent exponen- 
tial correlation length ^ becomes the appropriate reference length. 



Although the energy density (internal energy per unit volume) divided hj ks 
U{t,L,L)= -T^^i^^. (18) 

is completely determined by the free energy density f(t, L, L) it turns out 
that a separate discussion of the energy density is warranted because of its 
important role played in the mean spherical model in Sect. IV. From Eqs. 



([Tl| ) - (|T5D one obtains the prediction 

f/(t, L, L) = Us{t, L, L) + Uns{t, L, L) (19) 
where the singular part 

[/,(t, L, L) = T, i,''" U{L/i, L/L) (20) 

has the universal scaling function 

W(x, s) = -vx^-^l" dT[x, s)/dx (21) 

= -dv R+ x'^-i/'^ - 2{d - \)v x^-^-^l' - z/x^-i/^ dG{x, s)/dx (22) 

and where the leading nonsingular part 

Uns{t,L,L) = Uo{t) + 2U,{t)/L (23) 
has a regular t dependence with Uo(t) = —T'^dfo(t)/dT and 

U^{t) = ~T^d-^i{t)/dT. (24) 
For the surface energy density, Eq. (|^) implies asymptotically 

Usurfaceit) = 9 f sur f aceif) / dT (25) 
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= U,{t) - T, ^-'/'^ {d - l)v r^^— '^^/^ • (26) 

In Eqs. (^) and (^) we have used the hyperscahng relation 



dv = 2 — a 



(27) 



These scahng predictions have been confirmed by several field-theoretic re- 



normalization-group (RG) calculations of /^(t, L, oo) ||25| and of U. 



' surface 



5|, based on the ip^ continuum Hamiltonian with Dirichlet boundary 
conditions for the field v^(x). All calculations, however, were carried out 
within the dimensional regularization scheme which neglects cutoff effects. 
As pointed out by Dohm 35], an unresolved feature of the dimensionally 
regularized perturbative results for Usurface 0, 0, is a pole term 
~ (c? — 3)~^ that diverges in three dimensions. 



We note that the critical exponent of Ugurfaceif) 

1 — a — v={d— l)v — 1 



(28) 



is positive for ordinary critical points of the 0{n) universality class with 
d > 2 which implies a finite critical value Ugurfacei^) = Ui{0). By contrast, 
for the Gaussian model, 1 — a — u = (rf — 3)/2is positive only for d > 3, thus 
Usurfaceif) diverges for t ^ in (i < 3 dimensions (see Sect. III). 



We shall also consider the specific heat (divided by A;^ 

dU{t,LX) 



C{t,L,L) 



dT 



Cs{t,L,L) + Cns{t,L,L) 



From Eqs. ([19|) - (|23D we obtain the predictions 

C.(t,L,L) = i,"" L^'^C{L/i,L/L) 

and 



(29) 



(30) 



Cnsit, L,L) = dUoit)/dt + 2L-'dUi{t)/dt 



(31) 
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with the universal scahng function 

C{x,s) = V x^-^/" dU{x,s)/dx . (32) 

The scahng structure imphes that the surface specific heat Csurface{t) — 
dUsurface{t)/dT has a divergent singular part, 

C^urfaceit) = eo~'^S,3n./ace ^'"^ + dU^{t)/dT (33) 

with the surface scaling exponent 

as = a + 1/ (34) 
and with a universal amplitude 

Ai,surface = - {1 - a - u){d - 1) U At^^^^^ . (35) 

Finally we recall the prediction for the asymptotic scaling form of the sus- 
ceptibility 

xit,L,L) = Xb{t)fAL/^,L/L) (36) 

according to Eq. (1) where Xb(t) = x(^i oo) = ^''^ ^^'^ ^^^^ suscep- 
tibility. For L ^ ^, L ^ ^, the scaling function is expected to have the 
expansion 

/^(L/e, L/L) = 1 + c^C/L + 0{eiL\ e-^/^) (37) 
with the universal coefficient c^. For i > this implies 

X.nr/ace(t) =2im |[x(t,L,L) - Xbit)] ^| = A (38) 

with the surface scaling exponent 

75 = 7 1/ (39) 
and with the surface amplitude 



K,surface = ^ A ^0 ■ (40) 



14 



For T — > Ted the small L/^ behavior of the scahng function is expected to 
be 

f^{L/^,L/L)^B^{L/L){L/irl^ (41) 

with a finite universal amplitude B^{L/L) > which implies 

X(0, L, L) = ^-^^"B^iL/L) L^/^ . (42) 

In the following we examine the range of validity of these predictions for 
the exactly solvable Gaussian and mean spherical models in 2 < o? < 4 
dimensions. 
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III. Gaussian lattice model with free boundary 
conditions 



A. Lattice Hamiltonian 



We consider N continuous scalar variables (fj, ~oo < (fj < oo, on the lattice 
points Xj of a simple-cubic lattice with a lattice spacing a in a finite rectan- 
gular L X L'^~^ box of volume V = LL'^~^ = Na'^. We assume a Gaussian 
statistical weight ~ exp(— i/) with the lattice Hamiltonian 



H 



^0 2 



J 



E 2fi^ (^^ - 



<ij> 



(43) 



with a nearest-neighbor coupling J > 0. The factor (kBT) ^ is absorbed in 
H. The dimensionless partition function is 



n/^ 



^l-<i/2 



exp(-ii') . 



(44) 



In the bulk limit L — > oo, L — > oo, this model has a critical point at ro = 
for arbitrary d > 0. We assume that the temperature T enters only through 



T-Tr 



To = Oo 



, ao > . 



(45) 



A serious shortcoming of this model is the fact that it has no low temperature 
phase, i.e., no bulk limit exists for ro < 0. Nevertheless there exist nontrivial 
finite-size effects for tq > 0, as we shall see. 

We assume free boundary conditions in the d-th ("vertical") direction and 
periodic boundary conditions in the d—1 ("horizontal") directions. The d — 1 
"horizontal" coordinates and the "vertical" coordinate of the lattice points 
= (y ., Zj) are denoted by and Zj, respectively. The "bottom" and "top" 
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surfaces perpendicular to the vertical direction have the coordinates Zj = d 
and Zj = L, respectively, thus we have L/d layers of fluctuating variables. 
The variables in the bottom and top surfaces have only one neighboring layer. 
This is equivalent to assuming Dirichlet boundary conditions [ipj = 0) in the 
(fictitious) layers Zj = below the bottom surface and Zj = L + d above the 
top surface. The variables ipj can be represented as 

= Z-('^-i) (L + d)-^ J2 <^k,p exp(z k ■ y^.) V2 sm{p Zj) (46) 

with the Fourier amplitudes 

<^k,p = X] V'i exp(-ik ■ y^) sin{pzj) . (47) 



The sum 'Yli-^p runs over {d—1) dimensional k vectors with components fcj = 
2TT'mi/L, i = 1,..., d—1 with integers rrii = 0, ±1, ±2,..., in the range 
—7i/d < ki < n/d and over wave numbers p = Tcn/ {L + d), n = 1,2, L/d 
in the range < p < n/d. We see that, for L/d layers with free boundary 
conditions, the natural unit wave number in p space is 7r/(L + a) rather than 
71 /L. For each given p, there are {L/dY^^ variables (^k,p- Eq- ( ^6]) implies 
ifj = Q aX Zj = Q and (pj = Q ai Zj = L + d for arbitrary y^, thus we have 
a total number of = {L/d){L/dY~^ variables (^k,p- Substituting Eq. (^61) 
into Eq. (43) yields the diagonalized Hamiltonian 



2 

with 



H=l- L-^'^-^) (L + S)-i ('^o + + Jp) 0Kp 0~Kp (48) 



k,p 



4J 



Jk,d-i = — 5^ [1 - cos (kid)] , (49) 



i=l 



4J 

Jp = cos(pa)] . (50) 
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The Jacobian of the hnear transformation Lpj p of Eq. 



IS 



-N 



(51) 



Using Eqs. 
by ksT 



and (|5T1) we obtain the free energy density divided 



f{t,L,L) 



1 



2 



Invr 



(L/S 



In 2 



+ I L-^'-''> J2 In [(^0 + + Jp) ~a'] . (52) 

k.p 

In all calculations of this Section we shall keep the lattice spacing d finite. 

In the following we shall also consider film geometry (bulk limit in the d — 1 
horizontal directions). In Eq. (|52D, this corresponds to the replacement 

^"^''"'^Ek.p ^ EpL where = (27r)i-'^/ d''~'k with < n/d,t = 
1, 2, ...,d — 1, hence 

fit, L, cx)) = a-'^ Inn + ^L-' J2 f\n [(ro + J^-i + Jp)a^] ■ (53) 

P k 

A simplifying (but unrealistic) feature of the Gaussian model is that the 
critical point of the film system of finite thickness L is also determined by 
To = 0, i.e., it remains unshifted compared to the bulk critical point for all d. 
This differs from the case of the spherical model to be discussed in Sect. IV. 



B. Bulk properties 



First we briefly summarize some of the known bulk properties. The square 
of the second-moment bulk correlation length above is defined by 

e = lim lim i- S..(^->^.)'<^.^.> . (54, 

2a • < (pi (pj > 
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It is given by = Jq r^^ or 

e = eo = 1/2 (55) 

with 

^0 = {Jo/aoY^' , Jo = 2J . (56) 
From Eq. (^) we have the bulk free energy density for tq = a^t > 

1 1 f^^) 

f, = --d-Unn + -j In [(ro + J^,,) fi^] (57) 



q 



where /^''^ = (27r)-'^/ (/"^g with jg^l < 7r/d,i = l,...,d. Eqs. (0) - (^ 
are defined for all integer dimensions d = 1,2, ... They can be extended to 
continuous d, as usual, by means of analytic continuation via Euler's Gamma 
function F. From Eq. (|57D one obtains the singular part of fb for < d < 2 
and 2 < d <A 

hs = Rtr" (58) 
with the universal bulk amphtude 



F(3-ci/2) 
= 2^-Hnd-2) ■ 

The regular part of fb reads for < d < 2 and 2 < d < 4 

fo = a-'^ [ci + ro d^ C2 + d" C3 + ©(rjj d^)] (61) 

with (i-dependent constants Cj . The constants ci, £2 and C3 diverge for 
d ^ 0, d ^ 2 and c? — 4, respectively, where attains a logarithmic 
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dependence on roa^. The bulk susceptibility is simply Xb = = Jq 
which implies = a^^. The critical exponents are 



r] = 0, 



7 = 2z/ 



for all (i > 



(62) 



and 



a = (4 - d)/2 for < < 4 



(63) 



above Tc, in agreement with the hyperscaling relation Eq. ( pTf ) for < 4 
The prefactor in Eq. (^) is simply A^^C^q 



7/1/ 



The second-moment bulk correlation length ^ must be distinguished from the 
"exponential" bulk correlation length in the direction of the unit vector 
e = (xj — Xj)/|xi — Xjl which is defined via the large- distance behavior of 



anisotropic bulk correlation function G{xi — Xj) =< cpi ipj > For the 

special case where x = (x, 0, 0, ...) is directed along one of the cubic axes the 
correlation function decays exponentially as 1^ 



G(x) 



4J 



27[\x\ 



(d-l)/2 



sinh ( ^ 



{d-i)/2 



X e 



' [1 + 0(1 



(64) 



with the exponential correlation length 

^1 = — arcsinh \ — 



(65) 



We shall see that it is ^1 rather than ^ that determines the exponential part 
of the finite-size effects above T^. not only for periodic boundary conditions 
but also for free boundary conditions. 
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C. Free energy density 



In Appendix A we derive from Eq. (p2D the size dependent free energy density 
for box geometry for large L/a at fixed L/^ > and at fixed L/L for d > 1 



f{t,L,L) = h + 2fsurface{t) L-^ + Q{L/i,L/L)L-'' - hi-^L'-''\n2 



+ OihL 



where 



f surf ace ify 



dy[y-' [I + e-'y - 2e-'y h{2y)\ 







with the Bessel function of order zero 



TV 

1 f 

— d9 exp(2; cos( 
J 



Eq. (|66D contains the universal finite-size part 



Gix,s) = - dyy ^< {- 



sKis'y) 



d-l 



(66) 



(67) 



(68) 



1 /7r\(««-i)/2 
2\y) 



with 



K{z) = exp(— m^z) . 



(69) 



(70) 



We note that f surface depends on the lattice constant d, unlike the finite-size 
part Q{x,s). Using K(z) ~ (n/zY^'^ for z — > we obtain for film geometry 
(Z ^ oo) 



g{x,0) = ^ I dyy-' 



1/2 



'-K(y 

2 V4 



(d-l)/2 



(71) 
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The surface part remains of course identical with Eq. (|B7D. Eqs. (^) - (|7TD 
are apphcable to T = T,, and toT > T^. at fixed L/C,- The correct exponential 
large L/C, behavior of Q{L/C,, L/ L) at fixed T > Tc is not yet included in Eqs. 
( p9| ) and ( [7ID as it involves the "exponential" correlation length Eq. (|65D. 
For large L ^ at fixed T > Eq. ( |69D must be replaced by 

-id - l)(L/Z)('^+i)/2 [L/i2nC,)f e-^/«^ [l + Oi^l^'l-'^^)] . (72) 

Correspondingly, Eq. (p!6D must be modified for large L ^ ^. The nonuniver- 
sal last term ~ L^~'^ in Eq. (^) contributes to the regular part /„s(t, L, L) of 
/, thus Eq. ( p!5|) should be complemented accordingly. In order to clarify to 
what extent the d dependent term fsurfaceif) contains universal contributions 
we need to distinguish the cases l<(i<3,(i = 3, and 3 < c? < 5. For this 
purpose it will be useful to express the regular part linear in tq in terms of 
generalized Watson functions defined by |]T5l 



2tt 2-k 

WAz) = A-J-f '-^^ (73) 







oo 

J dye-^y[e-^yh{2y)Y . (74) 



1< d < 3 

For < ro ^ 1 and 1 < (i < 3 we obtain from Eq. (p7|) 

fsurfaceif) = fsurface{0) + ^turface ^ 

- b, ro d'-^J^' + 0(S/e") (75) 

with the universal amplitude 

r((3-rf)/2) 

2d+l T,{d-l)/2 _ 1) 



^surface nd+l „(d--lM2 / J i\ C^^) 
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and with the nonuniversal constant 



id-l 







This constant can be partially expressed in terms of generalized Watson 
function as 

bd = I [Wa-i{A)-2WM] 



+ ^Jdy[e-'yio{2y)Y 



oo 

1 1 {[e-^yioi2y)Y-' - (4vry)(i-^)/2} 



+ 



A 



(see Appendix A). We note that bfi does not depend on the arbitrary finite 
constant A > 0. 

The second term in Eq. (^) has the expected singular scaling form ~ 
The first term fsurface{0) and the third term ~ bd contribute to the regular 
part 2\I'i(t)L^^ of fns{t,L,L). Thus the surface contribution is in accord 
with the predicted universal scaling structure of Eqs. ([13|) , ([15|) , (|T^) and 
([I8|). We expect that bd depends on the lattice structure (see also Subsect. 
G). We note that both A^^j^ce ^^"^ coefficient bd diverge for ci — > 3 such 
that 



lim 



U _ A + 

" surface 



(79) 



has a finite limit b (see Appendix A) . The explicit expression for the constant 
b is given in Eqs. ( |HI| ) and ( P^ below. 
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d = 3 



For < ro ^ 1 we obtain from Eq. (^) aX d = 3 

f surfaced) = fsurfaceiO) - (IQti)-' H^/a) - 1 Tq J^^ + 0(5/^=^) . (80) 

The analytic expression for the nonuniversal constant h is 

A 

l2 



h = \jdy[l + e-'y - 2e-^y h{2y)\ e"^^ [Io{2y)f 



oo 

+ ^-jdy{[l + e-'y - 2e-2^/o(2y)] e"^^ [h{2y)f - {AT^yY'] 

A 

+ (327r)-i (1-Cij-lnA) (81) 

where C^; is Euler's constant. This constant can be partially expressed in 
terms of generalized Watson functions as 

h = \ [1^2(4) -21^3(0)] 

o 

A 

+ \jdy[e-'yh{2y)f 



oo 

\ I {[e-'yio{2y)Y - (Any)-' 



+ 



+ {327i)-\l-CE-\nA) . (82) 

Note that b does not depend on the arbitrary finite constant A > 0. We 
expect that b depends on the lattice structure (see also Subsect. G). 



While the first and third terms of Eq. (|80|) contribute to the nonuniversal 
regular part 2 '^i(t)L~^ of fnsit, L, L) the second (logarithmic) term is clearly 
a singular contribution to the free energy density. This implies that, for ci = 3, 
Eqs. (0) and (|16D must be replaced by 

fs{t,L,L,d) = L-'^{L/^,L/L,~a/0 (83) 
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where 



T{Lli,LlWali) = HL/i.L/L) - (STr)-^ {L/^f ln(e/S) (84) 

contains a logarithmic nonscahng term that depends on the nonuniversal 
lattice constant a. We conjecture, however, that the coefficient — (Svr)^^ is 
universal, i.e., independent of the lattice structure (see also Subsect. G). The 
scaling part reads 

T{L/i,L/L) = (127r)-i {L/if + g{L/^,L/L) (85) 



where Q{L/^, L/ L) is given by Eq. ([69| ) for box geometry and by Eq. (|7TD 
for film geometry. Correspondingly the universal scaling prediction of Eq. 



( [T^ ) must be replaced by 

fs{t,L,L,~a) = 4r'- [i8n)-'r' H^/a)]L-' 

+ L--'giL/^,L/L) + 0{L~'), (86) 

with the universal bulk amplitude = — (127r)^^ but without a universal 
surface amplitude. In Subsection D we shall see that this logarithmic behavior 
is related to the vanishing of the surface critical exponent of the Gaussian 
energy density at = 3. 



3 < d < 5 



For 0<a/^<^lwe obtain from Eq. (|67|) for 3 < d < 5 

fsurfaceii) ~ fsurface(fi^ Bd Tg Jq 

+ ^turface + 0^') (87) 

with the universal amplitude 

^surface ^{d-l)/2 - l){d - 3) ^ ^ 
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and with the nonuniversal constant 



oo 

Ba = \j dy[[l + e-'y - 2 e'^^ h{2y)\ [e'^^ my)]"''] > . (89) 







This constant can be expressed in terms of generahzed Watson functions as 
Bd=l [Wd^m + Wd-,{4) - 2W,{0)] . (90) 

o 

We note that Wd^i{0) exhibits a pole ~ (rf — 3)~^ near d = 3 where it can 
be represented as 



Wa-i{0) = 22-'^7r(i-'^)/2(rf-3)-M(3-'^)/2 + j dy[e-^yh{2y)\''-' 



oo 

+ jdy{ [e-'yh{2y)Y-' - (47ry)(i-'^)/^} . (91) 



The right-hand side of Eq. (91) is independent of the arbitrary constant 
A > 0. The first two terms of Eq. (^) contribute to the regular part 
2 \I'i(i(:)L~^ of fns(t,L) whereas the third term has the expected singular 
scaling form ~ ^^"'^^ Thus the surface contribution is in accord with the 
predicted universal scaling structure of Eqs. ([T3|), (p!^, (|T6|) and (0). We 
expect that depends on the lattice structure (see also Subsect. G). We 
note that Eqs. (|88|) and (|89D are the analytic continuations of Eqs. (|76D 
and (0), respectively, from d < 3 to d > 3 and that both Aj^j^^/ace 
coefficient B^ diverge for d 3 such that 



lim 

d^3+ L 



B^-A+ r 

" surface 



-- b (92) 
has a finite limit b, Eq. (^) (see Appendix A). 
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D. Energy density 



Eqs. (0) and (|52D yield the Gaussian energy density (divided hj ks) 

Uit, L,L) = - ^ E{ro, L, L, ~a) (93) 

with 

E{ro,L,L,a) = Z'-" L'^ Yl i^o + J^d-i + JpV' ■ (94) 

k,p 

In the following we must distinguish the cases T = T^. and T > T^. Using 
Eqs. ([T8|) and (|66|) - (|70| ) we obtain for T > and for large L/d and large 
^ /a at fixed L/^ > and fixed L/Lin2<(i<4 dimensions 

U{t,L,L) = Ubit) + 2Usurfaceit) L'' + T,^^^ £ {L / ^ , L /L) L^-'' + O {d L'-''){95) 

where Ub(t) = —T'^dfb/dT is the bulk part of U(t, L, L). Near Tc the surface 
energy density is given by 

Trd?-'^ 



Usurfaceif) — oti EsurfaceiXo 5, Jq ) (9^) 



with 



oo 

E^^rfaceiz) = ^ { [l + C-^^' - 26-^% (2^)] [6"^% (2y)] ""^"^^ } (97) 







where z = tq Jq^. Eq. (|97D can be expressed completely in terms of 
generalized Watson functions as 

Esurfaceiz) = Wd-liz) +W,_i(z + 4) - 2Wdiz). (98) 

The universal function £(x,s) = —\ x^^dQ{x, s)/dx of the finite-size part 
reads 

^^/* { (^)''^i|./WM/^(|)-il 

^ 

- if^j \e-y^y^-' (99) 
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for box geometry and 

oo 

1 



£{x,0) 



87r2 



dy 



1/2 



'-K(y 

2 V4 



^\ i'i-l)/2 

yj 



(100) 



for film geometry. We note that Egurface depends on the lattice constant a, 
unlike the finite-size part £{x,s). Eqs. (|95D - (|100D are not applicable to 
T = Tc for d < 3 where the functions Ugurfacei^) and £^(0,0) diverge. For 
3 < d < 4, Eqs. (||) - ( [T00|) are apphcable to both T = and T > % 
at fixed L/C, > 0. The correct exponential large-L behavior in terms of 
at fixed T > Tc is not yet included in Eq. ( |100|) . It can be derived from 
g{L/Ci,L/L),Eq. © . 



In order to see to what extent Egurface contains universal contributions we 
need to distinguish the cases 2 < d < 3,d = 3 and 3 < d < 4. 



2 < d < 3 



For < ro a <^ 1 and 2 < < 3 we obtain from Eqs. (^51) and 



U. 



surface 



it) 



-id-l)At^rfacee-' + 2d 



)~3-dl 



+o{d/c'-'xm 



•d anu ^surface < givcu by Eqs. ([77|) and (|76D . Eq. (|101|) 



where bd and A 

implies a divergent surface energy density ~ -vvith a universal surface 

amplitude (1 — '^)i^^^r/ace > ^ with the critical exponent 

(rf-3)/2<0, (102) 



a — u 



in accordance with the singular finite-size scaling part of Eq. (|26|) . Thus, for 
2 < d < 3, Us(t,L,L) satisfies the scaling prediction of Eqs. (jlUj) and ( PP]) 
with the critical exponent {1 — a)/^ = d — 2 and with the universal scaling 
function for x > 



U{x, s) 



duR+x'-' - 2(rf-l)z/A+,^_x'^-3 + 2z/^(x,s) (103) 



28 



where the universal bulk amplitude is given by Eq. (pHI). The function 
S{x,s) diverges as ~ x'^"^ for x — > 0. This divergence is cancelled by the 
surface term which implies the finite limit 



W(0,s) = lim W(x,s) = Sdis) 



where 



TT 



2„.Mrf-l 



2 l^Kis'y)] 



(104) 



(105) 



for box geometry and 



\ (<i-l)/2 
71 \^ " 



1/2 



K 



(106) 



for film geometry. By a separate calculation at T = Tc we find from Eq. (p3| 



f/(0,L,L) = f/,(0) + T,Co 



£Js)L'-'' + 2h,~a'~''L 



A 



surface' 



£d and hd diverge for d — > 3. 



(107) 



in agreement with the scaling parts, Eqs. (|105|) and (|106D. We note that 



d = 3 



For Q < tq / a? <^ 1 and for = 3 we obtain from Eqs. (|2 



and (H) 



U{t, L, L) = Ub{t) + [2Usurface{t) + T, ^o^^ ^(^v/^, L/L)J + O{aL~%108) 
with 

[/snr/ace(t) = ^T, ^q^^ [(STr)-^^ (^/S) + 26 - (STr)-^ + 0(S/o](109) 
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where b is given by Eq. (^Tj). Thus, as a special property of the Gaussian 
model, there exists a logarithmically divergent surface energy density with 
an explicit dependence on the lattice spacing d. This could have been an- 
ticipated on general grounds because ofl — a — z/ = for d = 3. (This is 
parallel to logarithmic terms for systems with periodic boundary conditions 
at the borderline dimension where the specific heat exponent a vanishes [^] . ) 
Thus, Eq. ( p6D is not applicable and the universal scaling prediction for the 
singular part Us{t, L, L), Eqs. (pO|) - (^2]), must be replaced by 

Us{t,L,L,~a) = T,^,^ L-'U{L/^,L/L,~a/0 (HO) 

where 

U{L/^,L/L,~a/0 = U{L/i,L/L) + (87r)-Mn(e/a) (111) 

with the scaling part 

W(x,s) = -duR+x + £{x,s) - (167r)"^ (112) 

We identify the nonsingular part as 

UUt,Lj) = Ut{t) - UUt) + T,^^%L-' + 0{dL-^). (113) 

The function S{x,s) diverges logarithmically for x ^ 0. This divergence is 
cancelled by Usurface{t) which implies that Eq. ( |108| ) has a finite limit for 
t at fixed L and L 

U{0,L,L) = Ub{0) + r,eo~' {-(Stt)-^ L-^ ln(a/L) + [6(s) + 2^] L"^} 
with the universal constant 

oo 



A 
A 

^ ^dyl - \ [sK{s'y)f [K{y/A) - 1] 





+ (167r)"MnA - (167r)"^ [1- C^; + 2 ln(27r)] . (115) 
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Note that b{s) is independent of the arbitrary constant A. We have confirmed 
the vahdity of Eqs. ( pT^ ) and ([TT5|) for = 3 by calculating U{0, L, L) 
directly from Eq. (193|) with tq = 0. Thus there exists a logarithmic nonscaling 
L dependence of the energy density at Tc, with an explicit dependence on 
a. As noted above we conjecture, however, that the coefficients (Stt)"^ and 
— (87r)~^ in Eqs. (|109|) and ( |114| ) are universal, i.e., independent of the lattice 
structure (see also Subsect. G). In Sect. IV we shall see that the logarithms 
in Eqs. (|109|) and ( |114| ) are the origin of the logarithms appearing in the 
three-dimensional mean spherical model with free boundary conditions. 



3 < d < 4 

For < ro ^ 1 and for 3 < < 4 we obtain from Eqs. (E^D and 



Usurfaceit) = UsurfaceiO) - T, ^,\d - e'"' + 0(5^^-^) (116) 

with the finite critical value 

UsurfaceiO) = T, S^'"'^ 5, > (117) 



where ^^^jace > are given by Eqs. (|88D and (|89D . The singular 

second term ~ yields a divergent slope ~ for t — > 0+, thus 

Usurface{t) has a uonuuiversal finite cusp at t = for 3 < c/ < 5, in con- 
trast to the case d < 3. As a consequence, only the temperature dependent 
contributions ~ ^^"'^ and ~ S{L/C^,L/L) to the energy density 

[/(t,L,L) = ^„,(0,L) - T.e^V-l)^^./-^'"'^"' 

+ T,^^'S{L/^,L/L)L^~^ + OidL'-'') (118) 

have a universal scaling form. The nonuniversal critical value ?7sur/ace(0), Eq. 
117|) , increases for d ^ 3. It enters the finite energy density at Tc 



Uns{0,L,L) = Ut{0) + 2Usurface{0)/L + OidL'-") (119) 
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which belongs to the nonuniversal nonsingular part of U{t,L,L) and which 
has a nonscahng L dependence ~ L^^. This L-dependence is nonneghgible. 
This will have significant consequences for the mean spherical model in c? > 3 
dimensions to be discussed in Sect. IV. 



E. Specific heat 



Eqs. (1^) and (|93D yield the specific heat (divided hj ks) 

Tao j^-ili-dJ2{ro + Jm-i + JpT' ■ (120) 



c 



From the first term of Eq. ( |120D we find full agreement with the finite-size 
scaling prediction, Eq. ([30|), in 2 < d < 4 dimensions with 2/z/ = 4, ajv = 
4 — (i. Specifically we find, for large Lja and ^ja oX fixed L/S,>Q and fixed 
L/L, the universal scaling function for a; > 



oo 

y 





for box geometry and 



e 



-yx'^/i-K'^ (121) 



C(x,0) = ^,jdyy[^-j e-^V4.^g 6-^4 (122) 







for film geometry. The evaluation of the second term of Eq. ( |120[ ) can 
be taken directly from subsection D for the energy density and yields only 
subleading corrections to the first term of Eq. ( |120| ). For T > Tc, Eqs. ( |121| ) 
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and ( |122| ) can be decomposed as 

1 



4 

^ 



C{x, s) = --d{d- 2) x'-' + X''-' 



d/2 



327r4 



2 Vl/ 



and 



Cix, 0) = --did- 2) i?+ x'^-^ + 



oo 

1 



327r 



dyyiTr/yy'^-'y' 



(124) 







where the bulk part (first term) contains the universal bulk quantity R^, 
Eq. (pTD, and where the surface part (second term) has the universal surface 
amplitude 

^hurface = ' r((5 - d)/2) (125) 

= \{d-l){^-d) (126) 

with Aj^j^^jace given by Eqs. (|76D or (^8|), in agreement with the predicted 
structure of the surface specific heat, Eqs. (|33| ) and (^). Eqs. ( |124| ) - ( |126| ) 
can be easily confirmed by calculating the derivative dU / dT from Eqs. ( p5D 

-(HOD. 



Eqs. ( |123| ) and ( |124D do not yet include the correct exponential part of the 



large-L behavior at fixed T > which involves the exponential correlation 
length ^1 [0] and which can be derived from Eq. (|72D . 
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F. Susceptibility 



The thermodynamic quantity of primary interest in the mean spherical model 
in Sect. IV will be the susceptibility. Important steps in the calculation of its 
finite-size properties can be performed already on the level of the Gaussian 
model. For box geometry the susceptibility is defined by 



xit.L.L) = <V:Vj> ■ (127) 

Substituting the representation Eq. (i3) into Eq. (|127]) we find 



Ea^ r / sin(na) 



;2 



X{t,m) = -^^^^^^[l-(-l)"] cot2(pa/2) <^o,p ^o,p> (129) 

with n = p [L + a)/7i = 1, 2, ... , L/a. From the Gaussian Hamiltonian Eq. 
(|48| ) we have < <^o,p <^o,p > = {L + a)(ro + Jp)~^ which leads to 

x(t,L,L) = — ^ J: |1 - (-1)1 . (130) 

L{L + a)^ ^ ro + Jp 

We note that this expression is independent of the dimension d and of L 
which is a special property of the Gaussian model. In App. B we evaluate 
Eq. (|130|) for large L/d. For large ^/a at fixed ratio L/C, > we find the 



scaling form 

x(t,L,L) = Xtf{L/0 = L-'/'^^L/O (131) 

with 7/z/ = 2 and the universal scaling function 

00 

f^{x) = ^1 dyil-e-y^'/^') [Kiy)~K{4y)], (132) 
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where K{z) is given by Eq. ([7D| ) and 

$(a;) = J^'x~'f^ix), (133) 

with X = L/C,- The leading terms of fxi^) for large x are 

f^{x) = 1 - 2x-^ + 0{x-'^) . (134) 

For X — s> we find lim^_^o 3^^^/x(2^) = 1/12, thus the leading L dependence at 
T = Tc is in accord with the scaling prediction Eq. (^2|), with A^^q'^^'^ = Jq^ 
and 



= ^ , (135) 



independent of the shape factor L/L. In the limit L ^ oo at fixed T > Tc, 
Eqs. ( |131| ) - ( |134| ) are valid only up to a nonuniversal exponential contribu- 
tion ~ e^^/^i in terms of L/^i rather than L/C, 0. 



G. Continuum approximation with Dirichlet 
boundary conditions 



As a first step towards the cp^ field theory with Dirichlet boundary conditions 
at finite cutoff A in a confined geometry we briefiy consider the continuum 
version of the Gaussian model. The comparison with the lattice version will 
serve to distinguish universal from nonuniversal contributions. 

The Gaussian continuum Hamiltonian for the scalar field v'(x) = fiy, z) 
reads 



Hfi. 



field 



(fx 



^0 2 , 1 /V7 ^2 



(136) 



This corresponds to the replacements 2J — ^ 1 and Jk,d-i k^, in 



Eqs. (H 



and (|57|). The wave numbers of the Fourier 
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components of the field v^(x) are assumed to have a finite sharp cutoff A. 
Although there exists no real system with a sharp cutoff the sharp-cutoff 
procedure is of significant conceptual relevance as it may signal important 
physical effects in real systems with subleading long-range interactions 0, 
0- 

As a consequence of the sharp-cutoff procedure, the bulk correlation function 
G'(x) =< {p{x.){p{0) > has an oscillatory power-law decay above Tc. For 
the anisotropic cutoff \qi\ < A ,i = l,...,d we find the anisotropic non- 
exponential large- distance behavior 

G(x) = 2^A^-2 ^ ^-2 ^-2)^1 -Q sin^Ax^ ^ ^ ^^Hx|/5) ^ (^37) 

i=i 

in contrast to the exponential decay of the lattice correlation function, Eq. 
(|6^). Thus the sharp cutoff induces long-range correlations, as expected 
0, H, For an isotropic sharp cutoff |q| < A, see Ref. 

The bulk free energy density above Tc is for 2 < d < 4 and for < d < 2 

fbjield = ^ + fojield (138) 

with the regular part 

f0,f^eld = A" [Cl + ro A-2 C2 + A"^ C3 + ...] . (139) 

The constants ci, C2 and C3 diverge for d ^ 0, d ^ 2 and c? — 4, respectively, 
where fbjieid attains a logarithmic dependence on ro A~^ . 

In the following we assume a L x 00 '^^^ film geometry, with Dirichlet boundary 
conditions (p{y, 0) = (p{y, L) = at the top and bottom surfaces. The 
continuum version of Eq. (p3D is 

ff,^^,(t,L) = A'^lnn + ^L-'J2 [ In [{ro + + p') A''] . (140) 

p i 
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The sum Yip runs over wave numbers p = nn/L, n = 1,2, ... in the range 
< p < A, and the components fc, are restricted to |A;j| < A, i = 1, 2, . . . , d—1. 
For large LA at fixed L/^ > we find for 2 < < 4 

ffield{t-,L) = fb 

,field 

+ g{L/0 L-" + O^A"-^ L-^) (141) 

where 



oo 



= ^ j dy y-' (1 - e-^) exp(-|/ ro A-^) (142) 



and 



oo 



/2(roA-^) = ^ I dy e"^ ^(y)^-^ exp(-y tqA-^) (143) 



with 



1 

S{y) = ^jd^ exp{-q'y) . (144) 



-1 



For the universal function Q{L/^) see Eq. ([71]) . 

In order to exhibit the universal and nonuniversal contributions to fsurface{t) 
we need to distinguish the cases 2<(i<3,(i = 3,3<(i<4. For 2 < < 3 
we find from Eq. (|142|) 

fsurfaceit) = fsurfaceit) ~^ ^surface ^ 

- hd ro A'^-^ + Oirl A'^"^) (145) 
with the nonuniversal constant 

oo 

h, = \jdy [(1 - e-y) Siy)"-' - i2n)'-' in/y)^''-'^'] . (146) 
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The universal amplitude A'^urface < is given by Eq. ([701) . For d = 3 we find 
a logarithmic nonscaling term similar to that of Eq. (|80|) , 



fsurfaceit) = LrfaceiO) - {IGir)-' HH) ' b + OiA''' T") (147) 

with the universal prefactor (IGtt)^^ and with the nonuniversal constant 

oo 

= ^ + 11 dy{{l-e-y) [Siyf - i4ny)-']} . (148) 







For 3 < c/ < 4 we find 



d-3 

+ Kurfacee-' + 0(A-^ ^ (149) 



fsurfaceit) fsurfaceit) Bd Tq A 

surface 



with the universal amplitude A'l^^^^^^ > 0, Eq. (|88D, and with the nonuni- 
versal constant 



oo 



Bd = \j dy SivY-^ (1 - e-y) > . (150) 



As expected, ffieidit, L) and fsurfaceit) contain the same universal parts as 
the corresponding functions of the lattice model. The nonuniversal constants 
fed, i) and Bd, however, differ from the corresponding constants 6^, h and Bd 
of the lattice model. For d — 3, 6^ and Bd are divergent. 

As an additional effect of the sharp cutoff, there exists the nonuniversal 
contribution to the free energy in Eq. ( |141|) . For d > 2 this term is 
nonnegligible compared to the universal scaling term L^'^ but it has a 
regular dependence on ro ~ t and therefore can be considered to belong 
to the nonsingular part fnsit,L) of the free energy density. Nevertheless it 
yields a leading nonuniversal contribution ~ A*^"^ to the Casimir force 
at Tc, similar to that for periodic boundary conditions discussed in Ref. M. 
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We briefly summarize the results for the energy density Ufieid{t, L) as derived 
from Eqs. ([Till) - O- For 2 < < 3, we obtain 



(151) 



The singular part is in full agreement with the flnite-size scaling structure. 
For d = 3, the energy density reads for large LA at fixed t > 



Ufieldit, L) — Ubjieldit) + 2 tJsurfaceif) + ^ ^ ^(-^/O 

+ C>(AL-2) 



L 



-1 



(152) 



with 



Usurface{t) = -T, [(Stt)-! ln(AO + 2h - {Su)-^ + 0{A-^ r')J • (153) 
In the limit t — * at fixed L we obtain 

[//,e/d(0, L) = f/,,/.e/d(0) + T, eo"' [(87r)-^L-i In(AL) + (6 + 26)^^^ 



with the universal constant 6, Eq. ( |115|) , and the nonuniversal constant 6, 
Eq. ( 148 ). The logarithmic nonscaling behavior in Eqs. (|153D and ( 154 ) is 



parallel to that in Eqs. ( 109 ) and ( |114|) of the lattice model in Section III D. 
The prefactors (Svr)"^ in Eqs. ( |153|) and ( |154|) are the same as in Eqs. (|109|) 
and (|114|) of the lattice model and are expected to be universal. 

For 3 < < 4, Eqs. ( p.l6| ) - ( p.l9| ) remain valid if a?~'^Bd is replaced by 
A'^^^Bd- Thus the surface energy density is 

Usurfaceit) = Usurface{^) " T, i,\d - 1) i^-' + 0{A-^ e'") (155) 

with a finite critical value 



U. 



surface 



(0) = Te ^0^2 A'^-s Bd > 



(156) 
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The singular part ~ ''is in agreement with finite-size scahng but is sub- 
leading compared to the regular part, as expected from the lattice model. 

For completeness we note that the specific heat and the susceptibility of 
the Gaussian continuum model with free boundary conditions are in full 
agreement with finite-size scaling for 2 < < 4, with the same scaling 
functions as those in Sects. Ill E and F for the lattice model, as expected. 



H. Dimensional regularization 



The method of dimensional regularization has been employed in all previous 
analytic calculations of finite-size and surface effects within the </7^ theory 
with Dirichlct boundary conditions. This method neglects cutoff and lattice 
effects. This is justified provided that the leading terms are universal. This 
is the case, however, only for d < d* where d* is a certain upper borderline 
dimension. In the present context of the Gaussian model with Dirichlet 
boundary conditions there exist the following upper borderline dimensions : 
d* — for the bulk free energy density ft,, d* = 1 for the surface free energy 
fsurface, d* = 2 foT the bulk energy density Ub, d* = 3 for the surface energy 
density Ugurface, d* = A for the bulk specific heat Cf, and bulk susceptibility 
Xb, and d* = 5 for the surface specific heat Cswface and surface susceptibility 
Xsurface- The mcthod of dimensional regularization correctly accounts for the 
leading universal parts only for d < d* (where cutoff and lattice effects are 
negligible corrections) and provides an analytic continuation to d > d*. It 
does not correctly describe, however, the cutoff and lattice dependent terms 
for d> d*. The upper borderline dimension d* — 3 for the Gaussian surface 
energy density will play an important role for the three-dimensional mean 
spherical model in Sect. IV. 

We begin with the analytic expression for the bulk free energy density of 
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the Gaussian model in d dimensions within the dimensional regularization 
scheme lET 



h,d^m{t) = - 2-^-1 TT^'^/^ Y{-d/2) ■ (157) 



According to Eqs. ( p^ - (pOD, fh,dimif) indeed agrees with the singular part 
fbs{t) of the bulk free energy in 2 < (i < 4 and Q < d <2 dimensions. The 
neglected terms are just those of the regular part /q. The latter is ultraviolet 
divergent for c? > dimensions according to Eqs. (pT| ) and ( |139| ). Near 
d = 2 the right-hand side of Eq. (|157|) has a pole ~ (rf — 2)~^ and therefore 
does not capture the logarithmic temperature dependence of fbs(t) in d = 2 
dimensions. 

Next we consider the size-dependent free energy density of the Gaussian 
model fdim{t, L) for film geometry within the dimensional regularization 
scheme. We find for general d 

hUt.L) = + ^At^rface^'-' L-' + QiL/OL-' (158) 

where ^^^/ace given by Eqs. (^) and (1711) . An alternative repre- 

sentation is given in Eq. (6.3) of Ref. Eq. (|1 58|) indeed agrees with the 
singular part fs{t, L) calculated in Sect. Ill for 2 < d < 3 and for 3 < < 4. 
For d = 3, however, the singular part depends explicitly on the cutoff or 
the lattice spacing according to Eqs. ( p.47| ) or (|86|). This is reflected in the 
dimensionally regularized result Eq. ( |158D only as a pole term ~ (c? — 3)~^ 
arising from ^^^/ace- Thus Eq. ( |158|) does not correctly describe the leading 
singular temperature dependence ~ tint of the surface free energy of the 
Gaussian continuum and lattice model in three dimensions according to Eqs. 
dH) and (ra. 



Finally we consider the dimensionally regularized result for the size-dependent 
energy density above Tc 

+ S{L/OL'-''] . (159) 
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For T ^ Tc this yields 

f/rf,„(0,L) = T.^^^S^L^-" (160) 

as confirmed by a direct calculation at T = Tc. We see that these expressions 
fail at (i = 3 where ^4^^^^^^ and £d do not exist because of pole terms ~ 
{d — 3)^^ as noted already by Dohm Eq. ( |159| ) does not capture the 
logarithmic divergence ~ Int of the surface energy density for T ^ at 
d = 3 according to Eqs. ( p.09| ) and ( p.53| ), and in Eq. ( p.60| ) the leading size 
dependence ~ In L of f/ (0, L) for L ^ oo at = 3 is lacking, compare 
Eqs. (|l3) and ( |T54D . 



Also for > 3, Eqs. ( |159| ) and ( |160| ) are not satisfactory since they contain 
only terms that are sub leading compared to the finite energy density at Tc, 
Eq. ( p.l9| ). The latter term that exhibits a nonscaling L dependence is 
missing in Eqs. ( |159| ) and ( |160| ). 



As far as the yj^ field theory is concerned, it is not clear at present whether 
these shortcomings are only an artifact restricted to the Gaussian approxi- 
mation (corresponding to one-loop order of the ip^ theory) or whether there 
exist further shortcomings at two-loop order. For this reason we do not 
consider universal finite-size scaling to be firmly established for nonperi- 
odic boundary conditions since the earlier field-theoretic results of Refs. 
p3| , 0, |25|, |2^, ^ are based on a perturbation approach using Gaus- 
sian propagators within the dimensional regularization scheme. On the other 
hand we note that the singular parts of both the specific heat and the sus- 
ceptibility are correctly described for the Gaussian model in 2 < c? < 4 
dimensions by means of dimensional regularization. 
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IV. Mean spherical model with free boundary 
conditions 



Again we consider N scalar spin variables , — oo < 5*^ < oo, on the lattice 
points Xj of a simple-cubic lattice with a lattice spacing a in a finite rectan- 
gular L X L'^~^ box of volume V = L L'^^^ = N d'^. We assume the statistical 
weight oc e~^^ with 



I <ij> i ) 



with a nearest-neighbor coupling J > 0. The "spherical field" fi{T, L, L, a) is 
determined as a function of /3 = (ksT)'^ and of L, L, a through the constraint 

d'^-^Y^ <Sf>=N = {L/~a){LraY~^ . (162) 



For a = 1, Eqs. ( |161| ) and (162) yield the standard formulation of the mean 



spherical model [|T^]. Keeping a as an independent nonuniversal parameter 
will facilitate the distinction between nonuniversal and universal contribu- 
tions. 

In the following we assume the same boundary conditions as for the Gaussian 
model of Sect. III. Thus the spin variables Sj can be represented as 

= Z"(^-i) (L + a)-i ^k,p exp(z k ■ y^.) ^2 sin(p Zj) , (163) 

and the diagonalized Hamiltonian reads 

H=]^ L-^"-^^ {L + S)-i 5^ (/i + J^d-i + Jp) S^,p ^-k,p (164) 

k,p 

with the shifted spherical field 

fl = 2fi~2Jodd-^ (165) 
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with Jo = 2 J. For Jk,d-i and Jp see Eqs. ( PU] ) and (0). The parameter 
/i(T, L, L, a) is determined imphcitly as a function of T, L, L and a through 
the constraint equation (|162|) which now reads 



Li^d^-i + Jm-1 + Jp)"' = (3 . (166) 

k,p 



The susceptibihty for T > T^^^ is 

xiT,L,L,~a) = Pj^J2<^^^i> (167) 



^L{lTW ~ ^"^^"^ cot2(pG/2) < So,p So,p > (168) 



p 

E ii - (-1)"! .,rrf/'' , im 



L(L + a) ^ MT,L,L,a)+J, 



with n = p{L + a)/TT = 1, 2, . . . , L/a which is parallel to Eqs. ( |127D - ( |130| ) 



for the Gaussian model. A significant difference, however, is the dependence 
of X on d, L, L and a through /i(T, L, L, a). 



A. Bulk properties 

First we recall some of the bulk properties. The bulk susceptibility at finite 
wave vector q above Tc^d is defined by 

X6(q) = Jim p J2 <S,S,> e-^^-^-^--^) . (170) 



The Gaussian structure of H, Eq. ( |164| ), implies |^ 

Xb(q)"' = /ife + Jq,rf (171) 
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where 



Hb = /i(T, (X), (X), a) = XbiO) ^ = (172) 
is determined implicitly by 

a"-^ (3-' J ih + Jn,d)~' = 1 • (173) 



q 



Eq. ( |173|) is the bulk version of the constraint equation ( |166| ). The square of 



the second-moment bulk correlation length is determined by the susceptibility 



according to |^| 



= XbiO) ^ [X6(q)'']q=0 = JoXb. (174) 
Setting Xb^ = fi'b = yields the bulk critical temperature Tc^d 

1 „ --J 2 /" T-l 



= /3c,. = a'-' / J^i . (175) 

I^B-l-c,d J 

q 

We note that Tc^d is nonzero for d > 2 and limrf^2+ Tc^d = 0. It is well known 
that the bulk critical behavior of the mean spherical model belongs to the 



universality class of the n- vector model in the large-n limit ||4^, thus the 



vanishing of ^ at c? = 2 is expected from the Mermin- Wagner theorem |5B 



In order to elucidate the role played by the borderline dimension d = 3 for 
the confined system we extend our analysis to continuous dimensions in the 
range 2 < d < 4. Eqs. ( |173D and ( [175| ) can be combined as 



iPcd -P)= Xb' j [J^4 iJ^,d + Xb')] ' ■ (176) 



This leads to the asymptotic critical behavior above T^^d for 2 < d < 4 

T-T, 



Xb = A^r\ e = eot-% t = —^, (177) 
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where 



7 = 2/(c^-2), u={d-2)-\ (178) 

and 

a = 2-du = {d-4:)/{d-2) . (179) 
We note that these critical exponents can be considered as the Fisher-renormahzed 



| 46| Gaussian critical exponents 



'yCauss ^Gauss /-■ 

7 = -7-1 ' ^ = ' 

^ ^Gauss ^Gauss 



and 



OiGauss /T m\ 

a = — , (181) 

-L (^Gauss 

as expected from the general theory of constrained systems with the 
Gaussian exponents of Sect. Ill 

iGauss = 1, l^Gauss = 1/2, acauss = (4 - (i)/2 . (182) 

The amplitudes are = Q/Jq and 

eo = fi[A,/(/5e,, Jo£)]'/('-') (183) 

with e = A — d and the geometrical factor A^, Eq. (|60D. The factor A^^q'^^'^ 
in Eqs. (|^) and ( ^2|) becomes simply A^^^^^" = Jq^. 

For completeness we note that at the lower critical dimension d = 2, the 
asymptotic behavior of Xb and ^ for T — Tc,2 = is exponential |^ and is 
derived from Eqs. (|173|) and ( p.74| ) as 



^ = c a exp i2npJo) , (184) 

= a^j-^ exp (47r/5Jo) (185) 
46 



with 



c = 0.03125 . (186) 

The vahdity of universal finite-size scahng to be derived in Subsection C 
below for 2 < c/ < 3 is expected to hold also for d = 2 above Tc^2 = in terms 
of the correlation length Eq. ( |184| ). 



B. Film critical temperature 



For film geometry {L oo) Eqs. ( |166| ) and (|169D are replaced by 



^-1 f(i2 + Jk,,-i + Jp)-' = (3 (187) 

p i 

and 

x(r.L.oo.a) = ^-/^ E II - (-1)1 ^.rf'^it 1 ■ 

L[L + a) ^ , L,oo,a) + Jp 

Unlike the box geometry, the film geometry introduces a considerable com- 
plication in that for d > 3 the film system of thickness L has its own 
sharp critical temperature Tc^d{L) > different from the critical tempera- 
ture Tc^d = ^c,d(oo) of the dimensional bulk system. As shown by Barber 
and Fisher [|l^, T^^^L) > Tc^d{oo) for d> A. Here we shall show that this is 
true also for 3 < c/ < 4. 

The condition for criticality of the film system is x{T-, L, oo, a)~^ = 0. This 
condition is satisfied at a critical value = fie < where jlc{L) is determined 
by the vanishing of the denominator of the lowest-mode {n = 1) term in the 
sumofEq. ([TSSD , 



~^c{L) = -^fl-cos-^) . (189) 
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We note that fic{L) is independent of d. The leading large L behavior is 

j2^{L) = - Jo7v\L + d)-"^ + 0[d'^{L + d)-^] . (190) 

According to Eq. ( |187| ), the corresponding critical value of l3c,d{L) = [kBTc^d{L)]~^ 
is then given by 

Pc,d (L) = ^f^- + + Jpy' ■ (191) 

P k 

Separating the lowest-mode (n = 1) term we obtain for general d 

~d~2 f ~d—2 ^1°" c 

n=2 •' 



The first integral in Eq. ( |192| ) is directly related to the critical temperature of 
a d— 1 dimensional bulk system [compare Eq. ( |175| )] and is infrared divergent 
for (i < 3, hence Pc,d{L) = oo or |]l5| 



Tc,d(L) = for < 3 (193) 

for any finite L, as expected from the Mermin- Wagner theorem We see 
that for the film system of finite thickness the dimension d = 3 plays the role 
of a lower critical dimension di = 3 up to which T^^L) vanishes. Thus, at 
finite temperature and in 2 < < 3 dimensions, there exists only one type 
of critical behavior for large L near the bulk critical temperature Tc^d > for 
the (i-dimensional film system of finite thickness. 

An analysis of Eq. ( |191| ) for d > 3 is presented in Appendix B. We find that 
Tc^L) is enhanced above Tc^oo) for ci > 3 for sufficiently large L ^ d. 
This enhancement is most naturally expressed in terms of the dimensionless 
parameter 

A/3 = Jo [f^cd (oo) - (3,,d (L)] . (194) 
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For large L ^ a the result is 

A/? = 4Ba d/L - Cd {d/lY'^ + 0(5^^/2 ^^gg^ 

with the nonuniversal amplitude Bd > 0, Eq. (p9D, and the universal ampli- 
tude 



Svr 



Ki-)-l 



TT 



(d-l)/2 



(196) 



with Cd > 0. Thus there are competing effects on Tc^^L) from the scaling 
term ~ L'^~'^ = V-/'^ and the nonscaling term ~ L^^ . The leading terms of 
the fractional shift are 

Tc,d{L) - Tc,rf(cx)) 



TrAoo) 



add/L - Cd{d/LY~^ + 0{d''l^ L~'''^) (197) 



with the positive amplitudes 

ad = ABd { I dy [e~'y Io{2y)]' 



(198) 



-1 



Cd = Cd{ dy [e-'yio{2y)Y 



(199) 



where Jq is given by Eq. (68). The amplitude can be expressed in terms 
of generalized Watson functions, see Eqs. ( PI1[ ) and (^). For o? = 4, 04 agrees 
with corresponding amplitude of Barber and Fisher |15| . 



We see that the positive shift of T^^diL) for > 3 is proportional to the same 
amplitude Bd > that determines the finite cusp of the Gaussian surface 



energy density, Eq. (|115| ). Thus the nonscaling Gaussian cusp and the 
nonscaling enhancement of Tc^d{L) for film geometry are closely connected. 
In the next Subsection we shall find that the Gaussian cusp is also responsible 
for nonscaling finite-size effects on the susceptibility for both box and film 
geometry in the mean spherical model for d > 3. 
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C. Constraint equation and susceptibility 



The crucial question is whether and for which dimension d the susceptibihty, 
Eq. ( |169| ), attains the universal scaling form of Eq. (|36| ) for large L,L,^ at 
fixed d. This requires to first analyze the size dependence of jl implied by Eq. 
( |166| ). Up to a constant factor, the left-hand side of Eq. ( |166| ) has the same 



form as the right-hand side of Eq. (0) for the Gaussian energy density, thus 



the constraint equation ( |166| ) can be rewritten as 

E{fi,L,L,d) = Pd^-'^ (200) 



where E{ro, L, L,d) is defined by Eq. (pi]). It is clear that any nonscaling 
L dependence of the Gaussian function E{fi, L, L, d) will cause a nonscal- 
ing form of the L dependence of fl which, through Eq. ( p.69| ), will in turn 
cause a corresponding nonscaling L dependence of the susceptibility. This 
mechanism explains the existence of a borderline dimension d = 3 between a 
scaling {d < 3) and nonscaling {d > 3) regime in the mean spherical model as 
a consequence of the size dependence of the energy density of the Gaussian 
model, for both box and film geometry. More specifically, we can anticipate 
nonscaling power laws for d > 3 arising from the nonscaling size dependence 
of the finite cusp of U surf ace ~ E surface accordiug to Eqs. (|116|) - ( |119| ). 



It turns out that the most natural parameter is not /i but rather the shifted 
parameter 

A/i = Jq-i [fl - fx,{L)] (201) 



where ^c{L) is given by Eq. (|189|) . We note that A/x > for box geometry 



for any finite L. In Appendix B we derive the following expression for the 
large L behavior of the susceptibility, Eq. ( |169|) , at fixed A/iL^ — vr^ > 



oo 
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where K{z) is defined by Eq. (|7D|). Eq. ( f^02| ) is valid for general d. 



In Appendix D we analyze Eq. ( pOOp for large L/a aX fixed shape factor 
s = (L + a)/L>0 near the bulk critical temperature. For small t = {T — 
Tc,d)/Tc,d > we find for 2 < d < 4 



+ 25/1-1 ^4 (A^) (203) 



where Esurface{z) is given by the Gaussian surface function, Eq. (p7\) - The 
universal finite-size part reads for box geometry 




d/2 



+ e 



y/4 



K 



d-l I ^_y^2/^^2 



(204) 



We note that the term K{y/A) — 1 = 2 exp(— 72^^/4) comes from the 
modes with the free boundary conditions {z direction) whereas the term 
sK{s'^y) = sJ2m=-oo^^T>{~s'^^'^y) comes from the modes with the periodic 
boundary conditions {d — 1 horizontal directions). For film geometry {s = 
0) the latter term is reduced to {n/yY^'^, and the universal finite-size part 
becomes 



£d{x, 0) 



167r2 



dy { - 



(d-l)/2 



2t'- 



1/2 



+ e 



?//4 



K 



(205) 



Eqs. (|203|) - ( |205| ) determine A/x implicitly near Tc ^ for large L as a function 
of t, L, L and d. In the absence of the a dependent term ~ Egurface, Eq. 



would yield a universal scaling form for A/iL^. According to Eq. ( |202| ), this 
would then imply a universal scaling form for the susceptibility. From the 
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analysis of the Gaussian model in Sect. Ill we know, however, that different 
scaling and nonscaling terms arise from Esurface depending on whether 2 < 
(i < 3, (i = 3, or 3 < d < 4. 



1. Finite-size scaling in 2 < d < 3 dimensions 

The asymptotic form of Esurface{^IJ'0?) for small A/ia^ reads for 2 < d < 3 

+ 86, + O ((A/xfi2)(^-2)/2) ^ (206) 



compare Eq. ( |101| ). Substituting Eq. ( p06| ) into Eq. ( |203| ) we see that the 



dependence on a is cancelled. This implies, for a given shape factor s = L/ L, 
that A/i has the universal scaling form 

A/i = L-^ Md{L/^,s) (207) 

where the scaling function A^rf(x, s) is determined implicitly by 

_^(.-2)/2 ^ ^ 2eA,'UM'J\ s) + 2e{d - l)A,' A^^^^^^Mt'^'" -(208) 

Substituting Eq. ( p07| ) into Eq. (|202|) confirms the finite-size scaling predic- 
tion, Eq. (PBD, with the bulk susceptibility Xh = PJo^^"^ and the universal 
finite-size scaling function for 2 < < 3 

oo 



At T = Tc^d this yields the power law Eq. p^) with 7/z/ = 2 and with the 
universal amplitude 



00 

^21 „, /^2l 
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The scaling results of Eqs. ( [^09D and ( |210| ) are applicable also to film geom- 
etry (s = 0) where /^(x, 0) and By,{0) are finite quantities for d < 3. We 
note, however, that both B^s) and B^O) diverge for d ^ 3. Specifically, 



for (i ^ 3 we find from Eq. ( pi Op for box geometry at fixed s > 

B^s) 271-^-^3 - d)-^ (211) 

whereas for film geometry 

B^O) ~ 27r-^ (3/2)2/(3-'^) . (212) 

The different types of divergencies are the consequence of a mode continuum 
for film geometry and signal different forms of violations of finite-size scaling 
aX d = 3 for box and film geometry, as will be confirmed in Subsection C2 
below. 

At fixed t > we find from Eqs. ( p07| ) - ( p09|) the leading large-L behavior 
for both box and film geometry 

= jV2^-i [1 + 2e{d -l){d- 2)-M, + 0{e/L')] , (213) 

and 

X = Xb{l-Md-l){d-2)-'A-/At^^^^^^ + 2]^/L + 0{e/L')} (214) 



where A~^urface < is given by Eq. ([76|). Eq. ( |214| ) yields the surface 



susceptibility, Eq. (|38|) , with the scaling surface exponent 

7, =7 + z/ = 3/{d-2) (215) 
and the surface amplitude 

Atsurface = ^0^ ^0 Md ' l){d - 2)-^A, + 2] . (216) 
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2. Violation of finite-size scaling in d = 3 
dimensions 



We recall that at c? = 3 there exists no sharp transition in box geometry 
and in film geometry of finite thickness other than the bulk transition for 
L — > oo at Tc,3 > 0, thus there is no compeUing reason to introduce a shifted 
reference temperature or to introduce a physical length scale other than the 
d = 3 bulk correlation length. At ci = 3 the Gaussian surface function reads 
[compare Eq. ( |109| )] 

Esurfacei^l^d^) = -(47r)-Mn (Afxd^) + 8\b- (167r)-^l + ©(A/i^/^^) .(217) 

Substituting Eq. into Eq. (|03D yields 



A/i^/^L = tL/^o + In (A/i^/^S) + 1/2 - 8nb 

+ 87r^3(A//^L,s) + ©(A/i^/^a) (218) 

with Ss{x,s) given by Eq. ( p04| ) . Here we have replaced shy s = L/L for 
large L/a ^ 1. Substituting Eq. ( |218| ) into Eq. ( p02| ) yields the suscep- 



tibility. As a consequence of the logarithmic term of Egurface in Eq. ( |217|) 
we now have a logarithmic dependence on d in Eq. (|2T8D that cannot be 
neglected. One expects that this causes only a logarithmic deviation from 
finite-size scaling. This will be confirmed for box geometry but not for film 
geometry where we shall find a power-law violation of finite-size scaling. The 
origin for this unexpected geometry effect at c? = 3 comes from the different 
small X behavior of the finite-size part £z{x-, s) for s > and for s = 0. 

2.1. Box geometry 

For fixed s > we find from Eq. ( p04| ) the small x behavior 



£^{x,s) = ]- s^x-' - -L inx + 0(1) . (219) 

2 OTT 
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Here the first term ~ is the contribution due to the k = 0,p = 7r/(L + 
a) mode which is the lowest mode of the discrete mode spectrum for box 
geometry. Eqs. (|218|) and ( pl9|) yield the leading L dependence at T = Tc,3 



1 + O 



ln(L/S) 



(220) 



A/i = 47rs^L-2 [ln(L/a)]"^ 

Substituting Eq. (|220|) into Eq. (|202|) leads to the susceptibility at T = Tc,3 
X = 2n-^jQ\s-^L'^ \n{L/a) + 0{L^) . (221) 
At fixed t > we find from Eq. (|218|) the leading large-L behavior 

;iV2 ^ |i _ + _ 2] ^/L + 0{M^/~a)nVL')} . (222) 

For the susceptibility, Eq. ( |202| ), this implies 



X = X6 1 1 + 



2 ln(e/5) + IQTTb - 6 C/L + 0{[\n{^/a)]^^yL^' 

which corresponds to the surface susceptibility at d = 3 

1 



~X.sur face (^) 



Jo' e 



21n(^/a) + 167r6- 6 



(223) 



(224) 



The logarithms in Eqs. ( p21| ) - (|224|) constitute logarithmic deviations from 
universal finite-size scaling, with an explicit dependence on the nonuniversal 
lattice constant a. Thus there exists no universal scaling form for box geom- 



etry at (i = 3 in the sense of Eq. (36). This is the consequence of the upper 
borderline dimension d* = ?> for the surface energy density of the Gaussian 
model. 

2.2. Film geometry 

A separate analysis is necessary for film geometry in c? = 3 dimensions since 
Eqs. ( |220| ) and ( p2 1| ) do not have finite limits for s — ^> 0. The susceptibility 
is again given by Eq. ( P02|) where now A/iL^ is determined implicitly by 



A/i^/^L = tL/io + ln(A/ii/2a) 



1/2 - 87r6 + 87r^3(A//^L, 0) (225) 
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with 

S{x, 0) 



dyy-U 1-2 



1/2 



y/i 



K 



- 1 



(226) 



We find from Eq. ( |226| ) tlie small x behavior 

3 



^3(X,0) 



Stt 



Inx + E'a + 0(a; In a;) 



(227) 



where 



+ 




+ e' 



I//4 



+ — dy y-^ 



-(n2-l)y/4 



(228) 



is independent of the arbitrary constant A > 0. The absence of a term ~ 
in Eq. ( |227| ) is a consequence of the fact that for film geometry there exists 
a mode continuum, without a discrete lowest mode. At the bulk critical 
temperature T = Tc^a Eqs. ( |225D and ( |227[ ) yield the constraint equation in 
the form 



A//2^ = ln(A//2~) ^ ^2 + 87r(E3 - 6) 

-31n(A//2L) + 0[A/iL2ln(A//L2)] . 



(229) 



We recall that the first logarithmic term in Eq. ( P29| ) is the signature of the 
(upper) borderline dimension ci* = 3 at which the critical exponent l—a—u of 
the Gaussian surface energy density vanishes whereas the second logarithmic 
term in Eq. ( |229| ) is the signature of the (lower) borderline dimension di = 3 
at which the film critical temperature Tcfi{L) vanishes, in accord with the 
Mermin - Wagner theorem pH]. Both logarithmic terms can be combined as 



ln(A//2~) _ 3 1j^^^^i/2^) _ ln(a/L) - ln(A/iL2 



(230) 
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and, after substituting into Eq. ( |229| ), yield the solution with a power-law 
(rather than logarithmic) L dependence 

A/i = A^hL-'' [1 + 0{a^/^L-^/^)] (231) 

with the nonuniversal amplitude 

A^ = exp|l/2 + 87r(E3-6)} . (232) 

Thus A/i has a nonscaling size dependence at T = Tc^a. Substituting Eq. 
into Eq. (|202|) leads to the susceptibility at Tc 3 



X = {JoA^~ar'L^ [I + 0{d'/'L-'/^)] (233) 

which constitutes a strong power-law violation of the scaling prediction x ~ 
L^, Eq. (P^D, in contrast to the logarithmic deviation in Eq. ( [221| ) for 
box geometry. The same violation persists in the critical region > 1. 
The unexpected difference between Eqs. ( p33|) and (|221|) results from the 



difference between the discrete mode spectrum for box geometry and the 
mode continuum for film geometry at the lower borderline dimension di = 3 
above which the film critical temperature becomes finite. We emphasize that 
the existence of this borderline dimension di = 3 which causes the second 



logarithmic term in the constraint equation ( p25| ) is not restricted to the 
spherical model. It remains to be seen whether similar geometry-dependent 
effects exist at d = 3 also in other models of 0{n) symmetric systems with 
n > 2 and with free (or Dirichlet) boundary conditions. 

At fixed T > Tc^^ for film geometry, we find the same leading large-L behavior 
as already given in Eqs. ( p22|) - ( p24| ) for box geometry, with a logarithmic 
(rather than power-law) violation of finite-size scaling. In summary, it is not 
possible to write x £1 universal finite-size scaling form, in the sense of Eqs. 
([^) and (^), in the region T > Tc,3 for film geometry at = 3. 

2.3. Comparison with Barber and Fisher 
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The susceptibility of the mean spherical model in film geometry with free 
boundary conditions was calculated by Barber and Fisher [|15| by a different 
mathematical technique. For d = 3 they introduced a " quasicritical tem- 
perature shift" and found that, for large n = L/d, there exists a scaling 
representation in terms of a scaled temperature variable nAK 



n 



Xbf 



- X(nAK) 



with a shifted inverse temperature deviation from the d 
temperature Tc,3 

J J 1 



(234) 
3 bulk critical 



AK 



C 



kBTc,3 

1 



inn 



In n + CBp/n , 



BF 



71n2 
IGtt 



The scaling function was represented parametrically via 



Xiz) 



(2/t/)tanh(y/2)] 
ln[(sinh y)/y] 



(235) 
(236) 



(237) 
(238) 



and was plotted for z > in Fig. 4 of Ref. ||T5|. For finite L/d ^ 1 and at 
AK = 0, Eq. ( |235| ) defines a "quasicritical" |^ temperature T(L) > Tc,3 
where 



J 



J 



knTiL) 



SnL \ d 



Cbf d/ L . 



(239) 



We note that all thermodynamic quantities are smooth functions of T near 
T(L) and that there exists no physical criterion for defining T{L). The analy- 
sis of the susceptibility in terms of AK in Ref. [|15| was restricted to AK > 
and did not include the region Tc^a <T< T{L), C,/L ^ 1 which is of principal 
interest for testing the scaling predictions of Eqs. (|1]) and (0). 

Our dimensionless susceptibility x/^^^ corresponds to Xbf- Our explicit result 
for x/^'^j however, is at variance with that of Barber and Fisher since Eqs. 
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and ( 226 ) cannot be reduced to the simple form of Eqs. ( ^34| ) - 
(|238|) . Our result differs from that of Ref. [jl5[ both in the region T^ s < T < 
T{L) and in the region T > T{L). A unique analytic comparison between 
x/a^ and xbf can be made in the region L/C, ^ 1 at fixed T above T{L). 
Substituting into xbf the large- 2; representation according to Eq. (9.17) of 

iH 

- 2{87iz)-^ ln(167re2) + 0{z-^ \nz) (240) 



Ref. 



X z) 



inz) 



and rewriting the resulting expression in terms of the bulk susceptibility x^bf 
and the asymptotic bulk correlation length [compare our Eqs. ( |177| ), ( |183| )] 

kBTr3 ,-1 



we find 



with 



Xbf 



bulk 

Xbf 



1 + 



a 
8^ 



J 



2 In 



C 



BF 



2 + - In 2 
2 



IGtt 



c 



BF 



i 

L 



O 



L2 



1 



W^2(4) 



W^3(0) 



(241) 



(242) 



(243) 



_4 ^ ' 2 

While the leading logarithmic term ~ [ln(,^/a)]^/L of Eq. ( |242|) agrees with 
ours in Eq. ( p23|) , the leading correction term ~ ^/L differs from ours. We 

does contain the last two terms of Cbf 
3) are missing in Cbf, Eq. 



note that our constant b, Eq. ([_ 
but the additional integral expressions in Eq. 



(243). Our integral expressions come from the finite part of Wd-i{0) for 
d —* 3, after subtracting the pole term ~ (rf — 3)~^, see Eq. (pi]). We 
believe that our result for b is correct since it has been obtained both by a 
calculation directly at d = 3 and by a calculation at d ^ 3, taking the limits 
d 3+ and d 3—. A further analytic comparison between our x/^^ 
and Xbf can be made directly at T = T(L) where xbf is simply given by 
= (nV2J)X(0) with X(0) = 1/12 according to Eq. (9.18) of Ref. flSl. 
Our result at T = T{L) depends on b and differs from the simple result for 
Xbf- Thus we doubt the correctness of the previous result for xbf for 
free boundary conditions at d = 3. 
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3. Violation of finite-size scaling in 3 < d < 4 
dimensions 



From Eqs. (0), (|116| ) and (|117|) we have the Gaussian surface function for 

d>3 



~2N{d-3)/2 



(244) 



where A^^^j^^^ > and are given by Eqs. ( [88| ) and (p9[). Substituting Eq. 
(HI) into Eq. (|20|) yields 



(245) 



We see that the finite value of EgurfaceiS^) ^ Bd > causes a nonnegligible 
term ~ {L/dY~^ in Eq. ( p45| ) that depends on the lattice spacing d. This 
will imply power-law violations of finite-size scaling for box geometry. 

At fixed s = L/L > 0, Eq. ( 245 ) yields the L dependence at T = T^^^ for 
3 < d < 4 



A/i = ^ s'^-^L^-'^ [1 + {d'^-^L^-'^)] . (246) 

For the susceptibility at T = Tc^d this implies 

X = 4Jo^d^-'^BdS^-''L'^''^[l +0 (S^-^L^-'^)] , (247) 
in contrast to the scaling prediction x ~ L^, Eq. (|4^). 
At fixed T > Tc^d we find from Eq. (|245|) the leading large-L behavior 

-1/2 ^ ji/2 ^-1 [1 _ ^^(^/^)^/^ + Oii^/dY-'e/L')] (248) 

60 



with the nonuniversal function for 3 < c? < 4 



m,(e/a) = 2eid-2)-'A,' [2i?,(^/S)^-^ - (rf - ,^_J • (249) 
Substituting Eqs. ( p48| ) and ( p49D into Eq. ( ^02] ) yields the susceptibihty 

X = Xb{l + [2m,i^/~a)-2]^/L + 0{{^/~af-''e / L')} . (250) 
The corresponding surface susceptibihty is for ^ ^ d and for 3 < < 4 

Xsurface = \ Jq' [2md{^ / d) ^ 2] ^ ~ ~ f'^/^'^-^) . (251) 

This is in contrast to the scahng prediction Xsurface ~ Xfe ^ ~ jf;-3/(c!-2)^ Eqs. 
(M) and (M). 



We see that the amphtudes of the leading nonscaling terms are proportional 
to Bd, both for the size dependence at Tc^d , Eq- ( |247|) , and for the tempera- 



ture dependence above T^^d , Eqs. ( p5CI| ) and (|251|) . Thus it is the cusp of the 
Gaussian surface energy density that is the origin of the nonscaling effects in 
the mean spherical model for 3 < < 4 rather than an enhanced transition 
temperature which does not exist for box geometry. 

For film geometry in 3 < d < 4 dimensions a new analysis of our solution 
would be necessary since there exists a sharp critical temperature Tc^d{L) > 
Tc,d{oo). Here we only note that the leading size dependence at fixed T > 
Tc^d{L) for film geometry is the same as given in Eqs. ( 248 ) - ( p^51| ) for 



box geometry. It is expected that a full description of the crossover from 
the L-dependent film critical behavior near T^^diL) to the d dimensional bulk 
critical behavior near Tc (^(oo) would involve two different correlation lengths. 
Our solution for x does provide the basis for such an analysis of a dimensional 
crossover which, however, is beyond the scope of our present paper. 
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Appendix A. Gaussian free energy density 



In this Appendix we derive the asymptotic form, Eq. (|66|), for the free energy 
density of the Gaussian lattice model for box geometry. We start from Eq. 
). Using the representation 



In ^ = y dy y-\e~y - e-'y^ 



(Al) 



we rewrite the L dependent part 



A/(t, L, L, 5) = fit, L, L, ~a) - Mt) + - ^-^1/^)'-^ In 2 
of the free energy density in the form 



(A2) 



(A3) 



with To = roa / (2 J) and 



^L/d,L/d,y)= S{oo,y) - S{L/d,y) SoiL/d^y) (A4) 



where 



S{L/a,y) = (a/L) ^exp [-2y{l - cos/c)] 



(A5) 



TV 

S{oo,y) = (27r)^^ J dxexp [— 2y(l — cosx) 



(A6) 



SD{,L/a,y) = (a/L) ^exp [-2y{l - cosp)] . 



(A7) 



The sums and run over dimensionless wave numbers in the range 
— vr < k = 2TTdm/L < it and < J9 = Tidn/{L + a) < vr with integers 
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m = 0, ±1, ±2, . . . and n = 1,2, . . . , L/d, as is appropriate for periodic and 
Dirichlet boundary conditions, respectively. In determining the large L/d 
and large L/d behavior of A/ at fixed finite ratio L/L it is important to 
distinguish the regimes < y < yo and y ^ yo with 

L + d , , , 

yo = A8 
a 

in the integral representation ( |A3|) . Accordingly we split 

A/ = i d-' (A/i + A/2) (A9) 

where 

yo 

A/i = jdy y-' e-'"y ^L/d, L/d, y) , (AlO) 



00 

A/2 = j dy y-^ e-'^^y ^{L/d, L/d, y) . (All) 
yo 

First we derive the leading L/d and L/d dependence of A/i. Since cosfc is a 
periodic function the sum S{L/d,y) satisfies the Poisson identity [^] 



00 } ^ 

S{L/d,y)= (27r)-W rffce™/^exp[-2|/(l -cosA;)] (A12) 

N=-oo 



^(00, y) + 2e-'y Yl nNL/d, y) (A13) 



N=l 

with 



TT 

F{M,y) = {2n)-^ J dk e'^^^ exp{2ycosk) = hii^y) (A14) 
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(A15) 



where 



Im{.z) 



- I dOe 

TT 



2 COS 9 



cos(M^) 



(A16) 



are the Bessel functions of integer order M (see, e.g., 9.6.19 of Ref. see 
also Eqs. (3.6) - (3.10) of Ref. 0). For y < yo at fixed L/L, the large L/a 
behavior of F{NL/d,y) for > 1 is F{NL/d,y) ~ 0(e~^/^), thus 

5(L/a, y) = 5(oo, y) + 0(6"^/^) . (A17) 

In order to determine the leading L/a dependence of the sum SD{L/d,y) for 
y < yo we first derive a representation of the one- dimensional integral 

b 



I{a,b) 



f\x)dx 



(A18) 



in terms of summations. The derivation is similar to that in Eqs. (A. 21) 
- (A. 30) of Ref. We assume the arbitrary real function f{x) of the 

real variable x to be well behaved in the interval a < x < 6, in particular 
we assume that f{x) has a convergent Taylor expansion around any x in 
this interval. We split the interval a < x < b into subintervals of length 
Ax = {b — a)/N > between the points Xi = a + iAx, i = 0,1, ■ ■ ■ , N, with 
Xq = a,XN = b. The integral I can be represented as 



I(a,b) = ^ f{x)dx . 

i=0 Z. 



(A19) 



For each interval we expand f{x) into a Taylor series around Xj+i (rather 
than around Xi as in Ref. |^ ) 



f{x)dx 



/(xi+i)Ax + 2^ 



^ (n + 1)! 

n=l ^ ' 



/("■)(x.+i)(Ax) 



rfa; (A20) 
(A21) 
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where /^"-'(x) = d"" f (x) / dx"^ . Thus we obtain 

/ f{x)dx = X;/(:^m)Ax + ^N^^^b) (A22) 



1=0 n=l 

where 

N-l 



K^;:\aM = $^/("n^m)Ax. (A23) 



1=0 



Since f{x) is an arbitrary function we may also apply Eq. ( [A22| ) to the 
function f'{x) instead of f{x). This yields an expression for K^^\a,h) m 
terms of higher derivatives, 

ir«(a,6) = f{h) - f{a) - E^^^^^ , (A24) 



n=l 



which can be substituted into the n = 1 term of Eq. ( |A22| ) . Successive 
application of this procedure permits one to express the difference 



I f{x)dx - 5^/(x,+i)Ax = RN{a,h) (A25) 

in terms of the differences of the derivatives at a and 6, 

^fik) ^ _ ^W(^) _ (A26) 

Note that RNia, b) differs from RN{ci, h) of Ref. |3ll. The result is 

MaM = -^[/(6)-/(a)] - ^^A/« + O ((Ax)^) . (A27) 

The coefficient of the O {{Ax)'^) term vanishes. Since As ~ 0{N~^) this 
representation is expected to converge rapidly for large if Af^'^^ remains 
sufficiently well-behaved for large k. Eq. (|A27|) differs from Eq. (A. 30) of 
Ref. by a minus sign in the first term. 
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We apply Eqs. ( |7^ - (|M7[) to the integral 

TT 

^(oo, y) = ^ j dp exp[-2y (1 - cosp)] (A28) 



where the integration variable p plays the role of x in the integral of Eq. 
( |A25| ). The sum on the left-hand side of Eq. ( [A25| ) now corresponds to 

1 ^ 

- ^ Ap exp I - 2y [1 - cos{nAp)] | = In^v) (A29) 

n=l 

with Ap = ir/N. Setting N = L/d we see that 

h/M = ^Sn{L/~a,y). (A30) 

We note that the derivative of the integrand of Eq. ( |A28|) with respect to p 
vanishes at p = and p = it. Eqs. (|A25|) - (|A30|) yield the leading large L/d 
behavior, for y < {L + d)/d, 

Sn{L/d,y) = S{oo,y) - ^ [l + e"^^ - 2e~'yio{2y)] + Oid^L') . (A31) 

In order to ensure that the discretization intervals Ap become sufficiently 
small for large L/d the restriction y < 0{L/d) was necesssary. For this 
reason, Eq. ( |A31| ) is applicable only to A/i, Eq. ( [A10| ), but not to A/2, Eq. 
([ATTl) . Using Eq. (|AT7D and substituting Eq. ( [A3l|) into Eqs. and 
(AlO) we arrive at 
yo 

^f^ = ^J dy {y-'e-'^y[S{^,y)t^[l + e-^y-2S{^,y)\] 


+ O (d^/L\ e-^/") . (A32) 
This can be combined with A/2 in the form 

A/1 + A/2 = ^jdy{y-'e-'^^y[S{^,y)Y-'[l + e-'y-2S{^,y)\] 


+ A/3(L/S,Z/a) + 0(a7L^ e-^/^) , (A33) 
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A/3 



dyy-^ e-'^^y\[S{^,y)Y{l + ~a/L) 



yo 



[S{L/d,yr-' Sn{Lra,y) - —[S{oo,y)] 



d-l 



(A34) 



The integral term in Eq. (|A33| ) represents the surface contribution of 0{L^^) 
to Af whereas A/3 will yield the finite-size part of 0{L~'^). Since y > yo 
in Eq. ( A34|) is sufficiently large it suffices to use the small k approximation 
—2y{l — cos A;) ~ —k'^y in Eq. (|A5|) , and similarly in Eqs. ( [A6|) and 



S{L/h,y)^{a/L)J2 ^'^'"^ = {a/ L)K{ATx^~a'L-^y) + ©(e""'^) , (A35) 



^(00,1/) ^TT-^ / dke-'^'y = {2n)-\n/yy/^ + ©(e""'^) 



SD{Lra,y)^Ca/L)J2 

p 

1 



(d/L) \K{7r'd\L + d)-'y) - 1 + ©(e"" ^) 



(A36) 



(A37) 



(A38) 



where K{y) is given by Eq. (^). Furthermore it is useful to turn to the 
integration variable 



z = An a y / {L + a 



(A39) 



Instead of yo we then have the lower integration limit Zq = An'^d/ (L + d) — 
for large L/d. This leads to 



~dr-i 



A/3 



(L + d 







1 /7r\('^-i)/2 
with the shape factor 



dz z 



exp 



-1 



d/2 1 



ro(L + a 
8J^ 



L + a 
L 
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[sK{§^z)]''-^[K{z/4) - 1] 



+ [1 + 0(sVL^)] (A40) 



(A41) 



For a we finally obtain Eqs. (p^ ) - (pUl). 



The surface free energy, Eq. (p?!), can be expressed in terms of the generalized 
Watson function, Eqs. (|73D and (0), as follows 



oo 

~l-d 



fsurfaceit) = ^ j dz [W,.,{z) + W,.,{z + 4) - 2W,{z)] . (A42) 



2 r-1 



It can be shown that for 7^ 3 there exists the following common represen- 
tation of the coefficients 6^, Eq. ([77|), and Bd, Eq. (p9|), of the regular term 
of f surf ace hucar in ro 

&d = 5, = ^ [l^d-i(4) - 2WM] 

o 

id-l 



+ i 1 rfy [e-^^/o(2^)] 


00 



"a 

+ 2-'^-V(i-'^)/2(d - 3)-i . (A43) 

This expression is independent of the arbitrary constant A > 0. For d —>■ 3+ 
and d —>■ 3—, the first two terms have a finite limit [W^2(4) — 2iy3(0)]/8 
whereas the last term exhibits a divergence ~ (d — 3)^^ that originates from 
Wrf_i(0)/8 for d 3+ according to Eq. (^TJ). The same divergence is 
contained in ^j"„r/ace' see Eqs. (^) and (|92D. 
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Appendix B. Susceptibility 



We rewrite Eq. ( |130|) as 
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X 



3(^$:{2-ii + (-iri} 



cot2(pS/2) 
'^o + Jp 



(Bl) 



2(1 



7,^ 



Y — 

JoL{L + a) ^ (1-cosp) [ro + 2(1 - cosp)] 
2a^ 1 + cos q 



1 + cosp 



JqL{L + a) 



COS g) [fo + 2(1 — cos g)] 



(B2) 



with fo = ^QCi'^ / Jq- The sums and run over dimensionless wave 
numbers p = nan/ {L + a) with integers n = 1,2, L/a and q = 27Tam/ (L + 
a) with integers m = l,2,...,L/ (2a) where we have assumed that L/a is an 
even integer. Using the decomposition 



1 + cos X 



1 



'1 — cosx) [ro + 2(1 — cosx)] 

1 



4 

+ — 



fo + 2(1 — cosx) 
1 



.2(1 — cosx) ro + 2(1 — cosx) 
and applying the representation Eq. ( [AID obtain 



2a 



7,5 



X 



JoL^{L + a) 



dy 



with 



^{L/d,y) = SD{L/h,y) 



i SD{L/a,y) 



where SD{L/a,y) is given by Eq. ( [A7| ) and 



2a 



SD{L/a,y) = — ^ exp[-2?/(l - cosg)] . 



(B3) 



^{L/d,y) (B4) 



(B5) 



(B6) 
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We distinguish the regions < y ^ yo = {L + d)/d and y > yo. The large 
L behavior of So{L/d,y) in the former region is given by Eq. ( |A31| ), the 
corresponding behavior of 5*/) is 



SDiL/d,y) = (l + d/L) S{oo,y) - nd/{L + d) + 0{ay L') (B7) 

which follows from the Poisson identity (see, e.g., Eq. (3.6) of Ref. [^). In 
the region y > 2/o "we may use the approximation (|A37|) and 

oo 

SD{L/d,y) = (2a/L) J]exp [-47r2aVV(^ + a)'] + 0(e-'^'^) (B8) 

n=l 

While the contributions of the region y ^ yo are important for the large L 
behavior of x at fixed T > Tc, the contributions of the region y ^ yo yield the 
finite-size scaling behavior of x in the critical region L/d ^ 1, ^/d ^ 1 at 
fixed ratio x = L/^ > 0, including the leading terms of the scaling function 
for large x, as given by Eqs. (|131|) - ( |135| ) for the Gaussian model. 

The derivation of x from Eq. (|169|) is parallel to that given above, except 
that fo is to be replaced by 

fi/Jo = Afi - 7Ty{L + dy + 0{d^L-^) , (B9) 

where Afi is defined by Eq. (pOl]) . At fixed M = Afi{L + a)^ > this leads 
to the large L behavior 



oo 



A(3{L + dy f J 1 - exp[-(M - n^)y/n'] 



Jqt^^L J I M-7r2 



exp [-(M - vr2)y/7r2] I [K{y) - K{Ay)] (BIO) 



fi2 



{L + dy 

where K{z) is given by Eq. (|7|). For L > a this yields Eq. (|202|) . 
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Appendix C. Film critical temperature for d > 

3 



In the following we consider Eq. (|191| ) for d > 3. Subtracting 

Moo) = h"-^ j J (Jm-1 + jpy' 

p k 

and using the representation 



(CI) 



dy e 



■zy 



(C2) 



for z > we obtain 



2J[/?,,,(oo) - = j dy^L/~a,y) 





(C3) 



^{x, y) = [S{oo, y)f ^ <; 5(00, y) - Soix, y) exp 



2y [1 — cos 



TT 



l + X 



(C4) 



where S{oo,y) and SoiL/d.y) are defined by Eqs. (^) and (^) [see also 
Eq. ( |A15| )]. It is important to distinguish the regimes < y < yo and y ^ yo 
with I/O given by Eq. (^). Accordingly we split 



yo 



dy^{x,y) = / dy<^{x,y) + / dy ^{x,y) = Ai + A2 . (C5) 



yo 



In Ai we use Eq. ( |A3 1| ) . A treatment similar to that in Eqs. (|A32|) - (|A39|) 
leads to 



L \L + a 



i-2 



1 + ^) + OCa''/^ L-"/^) (C6) 
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with the nonuniversal amphtude B^, Eq. and the universal amphtude 







(d-l)/2 



(C7) 



with Cd > 0- The first term ~ L ^ in Eq. ( |C6|) has a nonscahng L depen- 
dence whereas the second term ~ has the scahng L dependence ~ L^^'^. 
Eqs. (0) - (^) lead to Eq. (|19|). Rewriting 

l + e-^y-2e-^yio{2y) = 2e-2Mcosh(22/) - Io{2y)] (C8) 

and using (see 9.6.39 of Ref. PU[) 



cosh(2) - Io{z) = 2h{z) + 2h{z) 



> 



(C9) 



we see that Bd is positive and finite for d > 3, thus Tc^d{L) > Tc^dioo) for 
L ^ d. Using the representation 



2JpcAoo) 



dy [S{oo,y)Y 



(CIO) 



we obtain the fractional shift of the film critical temperature as given in Eqs. 
(|197|) - ( |199| ). We have verified that 04, Eq. (|198|) , agrees with the corre- 
sponding amplitude of Barber and Fisher |T5[ at d = A which was expressed 
in terms of the generalized Watson function, Eqs. ([73|) and ([f4D. The am- 
plitude ttd = 4:Bd/Wd{0), Eq. ( |198| ), diverges for d ^ 3. This divergence is 
cancelled by the next term of 0{L?'~'^) in Eqs. ( |C6D and ( |197| ). 
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Appendix D. Constraint equation 

We start from the constraint equation ( |166|) for box geometry where we 



decompose = Jq^/U + fic{L) and subtract /3c,d in the form of Eq. (|C1|) . 
Furthermore we add and subtract 

/3(A/i) = Ci'^-^y j (j^^^ ^ ^ j^yl pi) 

k p 

This yields 

/?e,<i - /3 = Ml + M2 , (D2) 
Mi = /3(A/i) - a^'^L'-'L-' J] ( JqA/x + /2, + J^-i + J^)"' , (D3) 

M2 = a'^-'A/i y y (JoA/i + + ^p)"' (^k,d-i + Jp)"' • (D4) 

k p 



Using the representation Eq. ( P2D we obtain 

00 

Jo Ml = j dy ^{L/d,L/d,y)exp{-Aiia^y) (D5) 



where ^{L/a, L/d,y) is given by Eq. (|A4|) . Again we spht the integral in 
Eq. (|D|) as = /o"' + = h + h with yo = (L + d)/d. For large L/d 
and L/d we find 

^ ^1"^^ + 2^(00, y)][5(oo,y)f-^exp(-A/ia2|/) 



+ 0(e-^/^S'^/2L-'^/2) , (D6) 
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a I L + a\ 1 



L \ d Svr^ 



K 



s Z^^ 



'7C\d/2 

-2{- 

.z . 



(D7) 



with zq = 47r^a/(L + a) and s = {L + a)/L where K{z) is given byEq. @ 
For large L/a we can let 2:0 ^ in Eq. ( p7| ). Evaluating the integral in Eq 



([D^ ) for small A/i yields for 2 < < 4 

JoM2 = £-iA,(A/ia^)('^-2)/2 



^2^ 



0(A/ia' 
Eqs. (OlD - O lead to 

+ 2£a ((A/x)^/' (L + a), s) [(L + d)/~af-'' 



(D8) 



(D9) 



where Esurface{z) and s) are given by Eqs. (|97| ) and (p04|) . Multiplying 
Eq. by ^^^j ^^gj^g 



we obtain Eq. ( p03| ) for L ^ a. 



(DIO) 
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